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最 负 态 名 的 程序 设计 竞赛 


A Google Code Jam ( GCJ) 


Google 公 司 举 办 的 一 年 一 度 的 程序 设计 竞赛 。 
O 影响 力 最 大 、 参 赛 面 最 广 的 程序 设计 竞赛 。 
O 每 道 题 都 有 Small 和 Large 等 不 同 规模 的 测试 数据 。 
© 自由 选择 喜欢 的 工具 在 本 地 运行 程序 ， 并 提交 结果 。 


A TopCoder 


TopCoder 人 公司 举办 的 程序 设计 竞赛 。 既 有 周期 举办 

的 SRM， 文 有 一 年 一 度 的 TCO。o 

O 和 在 75 人 分钟 内 措 战 难度 和 分数 递 增 的 三 道 题 。 

Oo 每 道 题 的 得 分 随 用 时 闻 减 。 

O 狸 有 的 挑战 险 段 ， 通 过 寻找 他 人 程序 的 漏洞 赚 取 宪 
引 的 分数。 

O 系统 评测 的 结果 和 在 最 后 才 人 公布。 


A ACM-ICPC 


美国 计算 机 协会 (ACM) 主办 的 面向 大 学 生 的 竞赛 -。 
O 历史 最 悠久 、 最 负 辟 名 的 程序 设计 竞赛 。 

O 二 名 选手 共用 一 人 台电 脑 的 团体 比赛 。 

O 5 个 小 HB 寺内 的 E 上 10 道 左右 的 问题 |。 


译 者 序 
3 


程序 设计 苋 赛 因 其 涉及 的 知识 面 广 ， 比 赛 形 式 激 烈 有 趣 ， 吸 引 了 越 来 
越 多 的 学 生 参 与 其 中 。 参 赛 


4 


者 不 但 可 以 从 中 锻炼 算法 设计 能 力 ， 还 能 够 提高 代码 编写 能 力 。 其 中 
WERE y PAGAS 


国际 知名 公司 的 重视 和 欢迎 。 
5 


本 书 的 几 位 作者 是 世界 公认 的 顶尖 选手 ， 在 竞赛 和 学 术 领 域 都 取得 了 
令 人 瞩目 的 成 束 。 他 们 结合 


自己 的 专业 知识 和 比赛 经 验 ， 将 自己 的 心得 和 技巧 集结 成 书 。 
6 


全 书 将 不 同 的 算法 和 例题 掖 专题 编排 成 小 下， 再 将 不 同 的 小 节 由 易 到 
难 分 成 四 草 ， 这 样 即便 是 初 


7 


出 茅 庐 的 新手 也 不 会 有 太 大 的 阅读 障碍 。 书 中 涵盖 了 在 程序 设计 苋 赛 
中 会 用 到 的 大 多 数 算法 和 技 


巧 ， 并 在 附录 中 补充 了 书 中 未 介绍 但 也 比较 有 用 的 算法 。 在 题材 的 安 
HE, FERES, EX 


分 明 ， 循 序 渐进 ， 不 以 华而不实 的 奇 技 淫 巧 误导 读者 ， 又 具有 一 定 深 
度 ， 相 信和 即便 是 经 验 丰 是 的 
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老将 同样 能 从 书 中 有 所 和 斩获。 本 书 在 结合 例题 进行 讲解 时 ， 不 是 简单 
地 堆砌 问题 和 代码 ， 而 是 注 


重 引 导读 考 更 好 地 理解 和 运用 算法 来 分 析 解 决 问题 。 对 于 正在 学 习 数 
据 结 构 与 算法 的 读者 而 言 ， 


9 
把 它 作 为 一 本 练习 和 拓展 的 参考 书 也 是 很 好 的 选择 。 


本 书 在 日 本 广 受 好 评 ， 还 先后 在 台湾 地 区 和 韩国 出 版 。 近年 来 程序 设 
计 竞 赛 在 亚洲 发 展 很 快 ， 在 


10 


中 国 大 陆 也 出 版 了 不 少 相 关 书 籍 ， 但 鲜 见 高 质量 的 佳作 。 所 以 ， 在 读 
SEBS, BAAR AES 


迫切 希望 中 国 大 陆 也 能 引进 这 样 的 好 书 。2012 年 初 ， 我 们 通过 作者 的 
推荐 了 解 到 了 本 书 第 二 版 的 


11 


出 版 ， 一 些 前 于 们 踊 距 翻译 计算 机 专业 书籍 的 经 历 也 或 舞 了 我 们 ， 让 
我 们 萌生 了 亲 上 日 翻译 此 书 的 


念头 并 联系 了 图 灵 教 育 。 非 常 幸运 的 是 ， 图 灵 教 育 也 正 考虑 引进 此 
书 ， 于 是 有 了 今天 呈现 在 各 位 


12 
ise TA BW AT Tal ASP SCH © 


在 翻译 上 ， 我 们 力求 做 到 既得 重 国内 选手 的 习惯 
的 表述 。 在 修正 原 书 中 的 一 
13 


， 又 符合 计算 机 专业 
些 笔 误 的 同时 ， 加 入 了 一 些 译 痢 注 ， 以 方便 国内 读者 理解 。 但 由 于 详 
者 水 平 有 限 ， 不 足 之 处 在 所 
14 


难免 ， 还 望 读者 多 多 包涵 ， 并 不 音 提出 意见 和 建议 。 
2 译 者 序 


们 的 一 些 疑 问 和 笔 误 给 予 了 


W TRR ERIAK], 


在 翻译 过 程 中 ， 秋 叶 拓 哉 、 知 田阳 一 和 北川 宜 秘 三 位 作者 耐心 地 对 我 


一 一 解答 和 确认 。 浙 江 大 学 的 陈 越 、 王 灿 和 俞 已 三 位 老师 不 但 将 我 们 
版 给 予 了 关切 和 文 持 。 在 此 
谨 对 他 们 表示 感谢 。 


还 拨 见 审阅 了 译 稿 并 提出 了 宝贵 的 意见 。 网 上 不 少 同好 也 对 本 书 的 出 
MFR ERI FREN 


2013 年 5 月 6 日 于 浙江 大 学 
IE 


如 今 ， 形 形 色色 的 程序 设计 竞赛 层出不穷 ， 听 说 过 Google Code Jam ` 
TopCoder、ACM-ICPC 的 读 


4 


者 八 介 不 在 少数 。 本 书 要 介绍 的 正 是 这 类 以 在 规定 时 间 内 、 又 快 又 准 
地 解决 尽 可 能 多 的 题目 为 目 


标的 程序 设计 竞赛 。 
5 


EAU E, HREZRZENETR, BATE LEPE E 
得 好 成 绩 也 绝 非 易 事 。 要 在 


程序 设计 竞赛 中 取胜 ， 不 仅 需 要 运用 灵活 的 想象 和 丰富 的 知识 得 出 正 
确 的 算法 ， 还 需要 一 气 呵 成 


地 实现 并 调试 通过 。 
6 


男 一 方面 ， 程 序 设 计 苋 赛 对 新 手 而 言 亦 非 吉 不 可 及 。 为 了 让 更 多 的 参 
赛 选手 体会 到 比赛 的 乐趣 ， 


大 多 数 比赛 都 会 准备 若干 面 加 初学 者 的 题目 。 另 外 ， 即 便 未 能 在 比赛 
中 取得 好 成 绩 ， 通 过 比赛 ， 
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也 能 够 使 目 己 的 能 力 得 到 有 效 的 锻炼 。 最 重要 的 是 ， 大 家 能 够 享受 到 
激烈 的 比赛 市 来 的 乐趣 o 


本 书 的 作者 们 参加 过 众多 程序 设计 竞赛 ， 在 平时 的 练习 和 学 习 中 ， 也 
获得 了 各 种 各 样 的 知识 与 技 
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巧 ， 本 书 将 这 些 知 识 技巧 总 结 成 有 册 ， 主 要 介绍 算法 及 其 在 相关 问题 中 
的 应 用 。 本 书 依照 由 易 及 难 


的 顺序 对 问题 进行 讲解 ， 章 节 的 编排 也 参考 了 主题 的 难 易 程度 及 其 相 
互 的 联系 ， 内 容 较 多 的 主题 


则 按 难 易 程 度 划 分 为 多 个 子 主题 分 别 介 绍 。 各 个 主题 由 算法 介绍 和 例 
古 讲 解 罕 插 而 成 。 


9 


只 要 是 具有 编程 基础 知识 的 读者 ， 均 适合 阅读 本 书 。 书 中 的 源 代码 均 
用 C++ 实 现 ， 不 过 只 用 到 了 


其 基本 功能 ， 所 以 即便 读者 不 熟悉 C++ 也 不 影响 阅读 。 
10 

【关于 再 版 】 
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令 人 惊喜 的 是 ， 本 书 的 第 1 版 受到 了 广大 读者 的 高 度 评 价 ， 在 此 表示 感 
W o Fl ERA 


PEP Fe AILS BAUR TAR > eA A RB MX AY AAA 
到 算法 ， 更 能 学 到 其 设计 和 


12 
运用 的 思想 。 这 正 是 本 书 划时代 的 腕 点 。 


本 书 第 2 版 退 加 了 计算 几何 、 搜 索 减 核 、 分 治 法 和 字符 串 相 关 算 法 4 个 
主题 。 此 外 还 退 加 了 方便 读者 


13 


加 深 理解 的 练习 题 ， 并 为 学 有 余力 的 读者 列 出 了 书 中 未 涉及 的 拓展 主 
题 ， 进 一 步 丰富 了 本 书 内 容 。 
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第 1 章 


鞭 势 待 发 一 一 准备 篇 


2 第 1 章 蕾 势 竺 发 一 -准备 篇 
1.1 何谓 程序 设计 竞赛 


首先 ， 让 我 们 来 说 明 一 下 程序 设计 竞赛 到 底 是 什么 。 


顾名思义 ， 程 序 设计 苋 赛 束 是 以 程序 设计 为 主题 举办 的 竞赛 。 世 界 上 
有 人 解 题 竞赛 、 性 能 竞赛 、 创 


意 苋 赛 等 各 种 各 样 的 程序 设计 苋 赛 。 本 书 主要 介绍 解 题 竞赛 。 


解 题 苋 赛 在 开始 时 会 告知 选手 题目 的 数量 ， 选 手 的 目标 是 解决 其 中 尽 
能 多 的 题目 。 > 程序 设计 竟 


赛 中 题目 的 形式 如 下 。 
抽签 


你 的 朋友 提议 玩 一 个 游戏 ， 将 写 有 数字 的 n 个 纸 片 放 入 口 激 中， 你 可 
以 从 口袋 中 抽取 4 次 纸 


片 ， 每 次 记 下 纸 片 上 的 数字 后 都 将 其 放 回 口袋 中 。 如 朱 这 4 个 数字 的 


Mim, Ma, E 


WU ERA o RAER THALE, ARMAR, TAR 
撕 破 口袋 ， 取 出 所 有 纸 


片 ， 检 查 上 自己 是 否 真 的 有 赢 的 可 能 性 。 请 你 编写 一 个 程序 ， 判 断 当 纸 
片上 所 写 的 数字 是 K 1, 


k2,...,kn 时 ， 是 否 存在 抽取 4 次 和 为 m 的 方案 。 如 果 存 在 ， 输 出 Yes 
; 人 否则， 输出 No > 


限制 条 件 


Ul<nz<50 


Sr: 


[] 1 <m< 108 
[] 1 < ki < 108 
FE 1 


m = 10 
k= {1, 3, 5} 

输出 

Yes (例如 4 次 抽取 的 结果 是 1、1、3、5， 和 就 是 10) 
1.1 何谓 程序 设计 竞赛 3 

样 例 2 

1 

输入 


n=3 


m=9 

k={1, 3,5} 

2 

输出 

No (不 存在 和 为 9 的 抽取 方案 ) 

3 

求解 这 个 问题 ， 可 以 编写 如 下 程序 。 
#include <cstdio> 


4 


const int MAX_N = 50; 

int main() { 

int n, m, KIMAX_N]; 

5 

/ 从 标准 输入 读 入 
scanf("%d %d", &n, &m); 

for (int i = 0; i < n; i++) { 
scanf("%d", &kli]); 

6 

} 

/是 否 找到 和 为 下 的 组 合 的 标记 
bool f = false; 

7 

/通过 四 重 循环 枚 举 所 有 方案 


for (int a = 0; a < n; a++) { 


for (int b = 0; b < n; b++) { 
for (int c = 0; c < n; c++) { 
for (int d = 0; d < n; d++) { 
8 


if (k[a] + k[b] + k[c] + k[d] == m) { 


f = true; 


} 

/ 输出 到 标准 输出 
if (f) puts(" Yes"); 
else puts("No"); 

10 

return 0; 

} 

11 


在 许多 比赛 中 ， 源 代码 一 经 提交 就 会 目 动 编译 并 运行 。 预 先 准备 好 的 
输入 文件 将 被 重 定 同 作 为 程 


4 第 1 章 蕾 势 竺 发 一 一 准备 篇 
Be 


当然 ， 程 序 的 运行 是 有 时 间 限 制 的。 在 大 多 数 比赛 中 ， 运 行 时 间 限 制 
在 若干 秒 。 一 旦 程序 运行 的 


时 间 超 过 了 限制 ， 程 序 就 会 被 强行 结束 ， 当 做 不 正确 的 解答 处 理 。 因 
此 ， 在 比赛 中 还 必须 考虑 高 


效 的 解法 。 


例如 ， 本 题 中 有 1 < n < 50 这 个 条 件 ， 像 上 面 那 样 单 纯 的 四 重 循环 的 程 
序 ， 不 用 1 秒 就 能 得 出 


Oy 
mt 


但 是 ， 如 果 变 成 1< n < 1000 又 会 怎样 呢 ? 四 重 循环 的 程序 即便 运行 很 
多 秒 也 不 会 结束 ， 这 将 被 


判 为 不 正确 。 不 过 ， 这 道 题 有 更 为 高 效 的 解法 ， 即 便 是 1 < n < 1000 的 
情况 ， 也 能 够 按 要 求 求解 


(将 在 1.6 节 中 再 讨论 ) 。 

由 此 ， 可 以 说 程序 设计 竞赛 是 综合 了 以 下 两 个 要 素 的 复合 竞赛 : 
口 设计 高 效 且 正确 的 算法 

0 正确 地 实现 
HB, 为 了 设计 算法 ， 

口 灵 活 的 想象 力 

O 算法 的 基础 知识 

也 是 必 不 可 少 的 。 

1.2 最 负 盛名 的 程序 设计 竞赛 5 
1 

1.2 最 负 盛 名 的 程序 设计 竞赛 


2 


o a TE. BUTDRST IR 
\ ne 


3 
1.2.1 世界 规模 的 大 赛 


它 是 Google 公 司 几 乎 每 年 都 会 举办 的 世界 规模 的 程序 设计 竞赛 ， 参 赛 
者 要 在 2~3 小 时 内 解决 大 约 4 


4 


道 题 。 一 旦 从 在 线 (Online) 进行 的 几 轮 预选 中 胜出 ， 就 能 够 参加 现场 
(Onsite) 总 决赛 。 该 赛 


事 的 特点 是 ， 每 道 题 都 备 有 Small 和 Large 两 组 输入 数据 。 即 便 是 难度 系 
数 较 大 的 问题 只 要 输入 


规模 足够 小 ， 依 然 可 以 简单 地 求解 ， 这 一 形式 深 受 广大 参赛 者 的 喜 
欢 。 另 外 ，GCJ 并 不 在 服务 器 


上 自动 执行 程序 ， 而 是 要 求 将 源 代 码 和 本 地 执行 的 结果 一 同 提交 > 
5 
1.2.2 回 高 排名 看 齐 


家 策划 并 举办 程序 设计 竞赛 的 公司 ， 它 举办 的 比赛 
涉及 多 个 领域 。 其 中 之 一 整 


6 


是 算法 (Algorithm) 比赛 ， 该 赛事 大 致 每 周 都 以 RM (Single Round 
Match) 的 形式 举办 一 场 ， 


HABI MRE ° 
(1) 
7 


Google Code Jam (GCJ) 


TopCoder 


在 1 小 时 15 分 钟 的 短 时 间 内 挑战 3 道 题 。 
(2) 提交 的 结 末 在 比赛 结束 前 古 不 知道 的 ， 整 个 过 程 中 稍 有 失误 ， 束 会 


AE AK O Zh e 


(3) 在 编码 阶段 (coding phase) 结束 后 ， 还 有 一 个 挑战 阶段 (challege 
phase) 。 该 阶段 可 以 查找 别 


人 代码 中 的 漏洞 。 如 果 能 够 提供 一 组 输入 数据 ， 使 别人 的 程序 返回 错 
RUTA, BACT EAN 


8 
外 的 分 数 。 


其 中 第 3 条 是 该 赛事 独一无二 的 特点 山 ， 也 是 阅读 别人 代码 的 好 机 会 。 
TopCoder A — (AZAR 


喜欢 的 等 级 分 系统 (rating system) ， 它 会 依据 SRM 的 结果 给 参赛 选手 
排名 。 男 外 ，TopCoder 还 会 


9 


举办 一 年 一 度 的 TCO (TopCoder Open) 公开 赛 。 一 旦 从 在 线 进 行 的 几 
轮 预选 中 胜出 ， 束 能 够 参 


加 在 拉 斯 维 加 斯 @ 举 办 的 总 决赛 。 


10 


O 随后 提 到 的 Codeforces 也 参考 TopCoder 提 供 了 类 似 但 不 完全 一 样 的 
hack 功 能 。 译 者 注 


O 最 初 儿 年 ，TCO 的 决赛 地 点 都 在 拉 斯 维 加 斯 ， 不 过 目 2011 年 起 ， 
年 的 决赛 都 选择 在 美国 不 同城 市 举办 ， 如 好 莱 坞 、 
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奥兰多 和 华盛顿 。 一 一 译 者 广 
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1.2.3 历史 最 悠久 的 竞赛 


ACM-ICPC 是 由 美国 计算 机 协会 (ACM) 主办 的 、 面 向 大 学 生 的 苋 
赛 ， 也 是 历史 最 悠久 的 程序 设 


计 竞 赛 。 这 是 一 个 三 人 一 队 的 团队 比赛 ， 选 手 要 在 5 个 小 时 内 解决 大 约 
10 道 题 。 因 为 比赛 中 三 名 


选手 共用 一 台电 脑 ， 题 量 又 比 其 他 赛事 多 ， 并 且 多 是 一 些 实现 复杂 的 
问题 ， 所 以 团队 配合 显得 异 


党 重要 。 想 要 从 日 本 参加 该 项 赛事 ， 首 先 要 参加 在 线 进行 的 国内 预选 
赛 ， 胜 出 后 才能 参加 亚洲 区 


域 赛 ， 取 得 前 几 名 的 好 成 绩 后 才能 够 参加 世界 总 决赛 。( 
1.2.4 面向 中 学 生 的 信息 学 奥 林 匹 元 驶 赛 


信息 学 奥林匹克 竞赛 是 学 科 奥 林 匹 克 竞 赛 的 一 种 ， 是 以 初中 生 和 高 中 
生 为 参赛 对 象 的 程序 设计 竞 


赛 。 在 日 本 ， 首 和 要 参加 日 本 信息 学 奥林匹克 竞赛 ， 取 得 优 寞 成 绩 
后 ， 才 能 作为 日 本 国家 队 选 手 


参加 国际 信息 学 奥林匹克 范 赛 。(2 其 他 比赛 都 需要 尽 可 能 快 地 解决 尽 


可 能 多 的 问题 ， 而 信息 学 奥 


林 匹 殉 竞 赛 只 要 在 规定 时 间 内 求解 问题 即 可 ， 成 绩 与 所 用 时 间 无 关 ， 
但 是 它 相 对 其 他 比赛 而 言 ， 


求解 每 道 题 所 伦 的 时 间 要 长 得 多 。 虽 然 是 面 癌 中 学 生 的 比赛 ， 每 年 所 
出 问题 的 难度 却 是 非常 高 的 。 


1.2.5 通过 网 络 自动 评测 


ACM-ICPC 


JOI-IOI 


Online Judge (OJ) 


在 互联 网 上 ， 有 一 些 被 称 为 Online Judge 的 系统 ， 它 们 能 够 自动 评测 以 
往 程 序 设计 苋 赛 中 的 题目 。 


利用 该 系统 束 可 以 练习 了 。 为 外 ， 其 中 一 些 Online Judge 也 会 定期 举办 
自己 的 比赛 ， 不 妨 去 参加 


一 下 。 在 此 列举 几 个 有 名 的 Online Judge ° 

[] PKU Online Judge (POJ) ——http://poj.org/ 

题库 中 有 大 量 的 题目 。 

O 会 津 大 学 Online Judge (AOJ) ——http://judge.u-aizu.ac.jp/onlinejudge/ 
还 包含 日 语 题 。 

[] Sphere Online Judge (SPOJ) ———http://www.spoj.pl/ 

允许 使 用 各 种 各 样 的 编程 语言 。 


[] SGU Online Contester 


http://acm.sgu.ru/ 
具有 模拟 参加 历史 比赛 的 虚拟 赛 功 能 。 


[] UVa Online Judge 


http://uva.onlinejudge.org/ 


老字号 Online Judge， 经 常 举 办 比赛 。 


[] Codecorces http://codeforces.com/ 


2 TopCoder 一 样 定期 举办 比赛 ， 又 同 其 他 网 站 一 样 不 断 维护 历届 题 
车 o 


O PRAM ASE RENEA REE, RS AE 
赛区 的 网 络 预赛 和 现场 区 域 赛 并 获得 前 几 名 。 当 然 


根据 规则 也 有 可 能 从 亚洲 其 他 地 区 获得 出 线 权 ， 其 具体 规则 比较 复杂 
并 可 能 不 断 变化 ， 大 家 可 以 从 网 上 获得 最 新 的 规 


则 。 一 一 译 者 注 


O 中 国 大 陆 的 中 学 生 首先 要 闻 过 全 国联 赛 (NOIP) 、 全 国 竞赛 
(NOI) MEZZNER (CTSC) 三 关 ， 才 能 参加 国际 信 


息 学 奥林匹克 竞赛 (ON) 。 一 一 译 者 注 

1.3 本 书 的 使 用 方法 7 

1 

1.3 本 书 的 使 用 方法 

2 

在 此 ， 束 本 书 所 涉及 的 内 容 、 使 用 方法 及 注意 点 做 一 下 说 明 。 
3 

1.3.1 本 书 所 涉及 的 内 容 


本 书 主要 讲解 程序 设计 竞赛 中 的 经 典 问题 和 基础 算法 ， 并 介绍 便捷 的 
实用 技巧 。 如 采 仅 仅 是 死记 


4 


经 典 问题 和 基础 算法 ， 遇 到 难 解 的 应 用 问题 或 是 需要 灵活 想象 力 的 问 
题 时 ， 仍 然 会 难以 下 手 。 因 


此 ， 为 了 加 深 理 解 ， 我 们 通过 选 目 POJ 的 经 典 题 和 部 分 原创 题 来 介绍 实 
EPA e 


男 外 ， 每 章 末尾 都 备 有 挑战 GCJ 中 实战 题目 的 小 栏目 ， 里 面 都 症 精 选 出 
来 的 题目 。 尽 管 要 找到 正 


5 


确 的 解法 恕 人 不 太 容 易 ， 还 是 建议 读者 匈 目 己 试 着 多 思考 一 下 。 在 此 
基础 上 再 阅读 题解 ， 能 够 得 


到 更 深刻 的 理解 。 


当然 ， 在 本 书 所 介绍 的 解法 之 外 ， 还 会 有 更 向 活 或 更 高 效 的 解法 。 大 
家 不 妨 多 试 独 去 思考 一 下 别 


6 
的 解法 e 

1.3.2 所 用 的 编程 语言 
7 


比赛 中 可 用 的 编程 语言 各 色 各 异 ， 而 C++ 在 几乎 所 有 比赛 中 都 可 用 。 它 
的 运行 速度 快 ， 库 函数 丰 


富 ， 因 而 人 气 很 高 。 本 书 选择 C++ 作为 所 用 的 编程 语言 ， 并 基本 按照 
g++ 的 规范 来 编写 源 代码 。 


8 
1.3.3 题目 摘 述 的 处 理 


在 世界 规模 的 大 赛 中 ， 理 所 当然 是 用 英语 来 描述 题目 的 。 不 过 ， 因 为 
题目 描述 中 的 英语 不 那么 难 ， 


所 用 的 单词 往往 也 非常 有 限 ， 所 以 很 快 束 能 习惯 。 当 然 ， 这 不 是 英语 
ei, FE te RF A A te 
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用 的 。 另 外 ， 其 中 有 些 比 赛会 针对 日 本 选手 提供 日 语 版 的 题目 描述 。 
英语 的 阅读 理解 不 是 题目 的 


关键 ， 因 此 本 书 的 题目 都 与 最 开始 的 例子 一 样 用 中 文 概述 4)。 
1.3.4 程序 结构 


10 


在 许多 比赛 中 ， 程 序 都 从 标准 输入 按 指定 格式 读 入 数据 。 输 入 并 非 问 
题 的 关键 ， 所 以 本 书 的 程序 


11 
© 原 书 为 用 日 语 搬 述 。 FAE 
871% ERRAR 


都 假设 输入 数据 已 经 由 main 函 数 谈 入 并 保存 在 全 局 变量 中 ， 再 通过 调 
用 solve 画 数 来 求解 。 例 如 


对 于 最 初 的 例子 ， 程 序 将 变 成 这 样 。 
/ 读 入 输入 数据 后 保存 在 这 里 


int n, m, kIMAX_N]; 


void solve() { 

bool f = false; 

for (int a = 0; a < n; a++) { 

for (int b = 0; b < n; b++) { 

for (int c = 0; c < n; c++) { 

for (int d = 0; d < n; d++) { 

if (k[a] + k[b] + k[c] + k[d] == m) { 
f = true; 


} 


} 

} 

} 

if (f) puts(" Yes"); 
else puts("No"); 
} 

1.3.5 练习 题 


每 章 末尾 都 会 介绍 与 本 章 所 涉及 主题 相关 的 题目 。 请 利用 它们 来 加 深 
理解 、 巩 固 知 识 、 培 养 实践 


能 力 。 各 个 主题 下 的 题目 大 致 是 按照 难 易 程度 排列 的 ， 其 中 亦 包含 非 
党 难 的 应 用 问题 。 


1.3.6 读 透 本 书后 更 上 一 层 楼 的 练习 方法 


独自 练习 提高 时 ， 不 妨 同时 使 用 Online Judge 和 TopCoder 的 Practice 
Room。 特 别 是 TopCoder， 既 


提供 了 解 题 教程 ， 又 可 以 阅读 别人 的 代码 ， 当 你 无 论 如 何 都 想不到 解 
法 时 ， 还 可 以 把 它们 作为 参 


考 ， 因 而 在 此 推荐 。 
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1.6 轻松 热 号 


本 下 通过 对 儿 个 问题 的 解析 ， 和 大 家 一 起 了 解 程序 设计 壳 赛 中 题目 的 
风格 ， 学 习 设计 算法 并 佑 


算 复杂 度 的 过 程 。 其 中 有 一 些 问题 比较 复杂 ， 不 能 想到 它们 的 解法 也 
不 要 紧 ， 能 够 通过 阅读 题 


解体 会 到 其 中 的 趣味 就 好 。 
1.6.1 先 从 简单 题 开 始 
三 角形 


A nthe, 棍子 i 的 长 度 为 ai。 想 要 从 中 选 出 3 根 棍子 组 成 周 长 尽 
可 能 长 的 三 角形 。 请 输 


出 最 大 的 周 长 ， 厦 无 法 组 成 三 角形 则 输出 0。 
用 5 根 棍子 组 成 三 角形 的 例子 
限制 条 件 


[]}3<n<100 


[]}1<ai< 106 
样 例 1 
输入 


的。 NETA CAMARO 
nan A 
= 8 
11T 


y 


AN 


1.6 轻松 热身 17 
输出 

1 

12 (选择 3、4、5 时 ) 
样 例 2 

2 

输入 

n=4 

a = (4, 5, 10, 20} 
3 

输出 


0 〈 无 论 怎么 选 都 无 法 组 成 三 角形 ) 


选择 3 根 棍子 ， 它 们 能 组 成 三 角形 的 充 要 条 件 为 


4 
最 长 棍子 的 长 度 < 其 余 两 根 棍 子 的 长 度 之 和 o 
5 
6 
7 


能 够 组 成 三 角形 的 条 件 


于 是 我 们 可 以 试想 这 样 一 种 算法 : 首 移 用 三 重 循环 枚 举 所 有 的 棍子 选 
FJR, BAH ERAH BE 


8 
否 组 成 三 角形 。 如 果 可 以 ， 那 么 该 三 角形 的 周 长 就 是 备 选 答案 。 


这 里 用 了 三 种 循环 ， 所 以 复杂 度 是 O(n3)。 将 n=100 代 入 n 3 得 到 
106， 可 知 这 个 复杂 度 是 足够 低 


的 Q)。 
9 
/ 输入 


int n, alMAX_N]; 
void solve() { 

10 

int ans = 0; // 答案 


/让 i<j<k， 这 样 棍 子 就 不 会 被 重复 选中 了 


for (int i = 0; i < n; i++) { 


O 本 题 还 有 O (nlogm ) 时 间 更 高 效 的 算法 ， 留 给 有 兴趣 的 读者 思考 。 


RESNE 
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for (int j =i+1;j <n; j++) { 

for (int k =j + 1; k < n; k++) { 

int len = a[i] + a[j] + alk]; // WK 

int ma = max(ali], max(alj], alk])); / 最 长 棍子 的 长 度 
int rest = len - ma; / 其 余 两 根 棍子 的 长 度 之 和 
if (ma < rest) { 

/ 可 以 组 成 三 角形 ， 如 果 可 以 更 新 答案 则 更 新 
ans = max(ans, len); 

} 

} 

} 

} 

/ 输出 

printf("%d\n", ans); 

} 


1.6.2 POJ 的 题目 Ants 
Ants (POJ No.1852) 


n 只 蚂蚁 以 每 秒 1cm 的 速度 在 长 为 工 cm WEF EMG AT o SHS eE 
子 的 端点 时 融会 摊 落 。 


由 于 竿 子 太 细 ， 两 只 蚂蚁 相遇 时 ， 它 们 不 能 交错 通过 ， 只 能 各 目 反 辐 
EHER o X TEREE , 


我 们 知道 它 距离 特 子 左 端 的 距离 xi ， 但 不 知道 它 当前 的 朝向 。 请 计算 
所 有 蚂蚁 落下 竺 子 所 需 


的 最 短 时 间 和 最 长 时 间 。 
笔 子 和 蚂蚁 的 情况 
限制 条 件 


[]1<L<106 


[]1<n<106 


[O<xi<L 
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样 例 

1 

输入 


L=10 


n=3 


x= {2,6,7} 


输出 


min=4 (`` A) 
max=8 (` `A) 


w 


IS — SAR ROBIE, IES BETRETEN 
组 合 ， 这 可 以 利用 递归 函数 


实现 (02.17) 


4 
每 n 只 蚂蚁 就 是 2x2x...x2=27m 种 。 如 
R n 比较 小 ， 这 个 算法 


是 可 行 的 ， 但 指数 函数 随 着 的 增长 会 急剧 增长 。 


gı 


2 n 增长 的 趋势 


10 20 

30 

100 10000 

1000000 

2n 232 1024 1048576 


109 1030 103010 10301030 


6 


穷竭 搜索 的 运行 时 间 也 随 之 急剧 增长 。 一 般 把 指数 阶 的 运行 时 间 叫 做 
指数 时 间 。 指 数 时 间 的 算法 


无 法 处 理 稍 大 规模 的 输入 。 


接 下 来 ， 让 我 们 来 考虑 比 穷竭 搜索 更 高 效 的 算法 。 首 先 对 于 最 短 时 
间 ， 看 起 来 所 有 蚂蚁 都 朝向 较 


7 


近 的 只 点 走 会 比较 好 。 事 实 上 ， 这 种 情况 下 不 会 发 生 两 只 蚂蚁 相遇 的 
情况 ， 而 且 也 不 可 能 在 比 此 


更 短 的 时 间 内 走 到 年 子 的 端点 。 


按 下 来 ， 为 了 思 芳 最 长 时 间 的 全 况 ， 让 我 们 看 看 蚂 晓 相 韶 时 会 发 生 什 
A o 


8 
9 
相遇 后 会 发 生 什么 
10 


事实 上 ， 可 以 知道 两 只 蚂蚁 相遇 后 ， 当 它们 保持 原样 交错 而 过 继续 前 
进 也 不 会 有 任何 问题 。 这 样 


11 
@ 也 叫 蛮 力 搜索 ,口语 中 常 简称 暴 搜 。 一 一 译 者 注 
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看 来 ， 可 以 认为 每 只 蚂蚁 都 是 独立 运动 的 ， 所 以 要 求 最 长 时 间 ， 只 
求 蚂蚁 到 和 合子 端点 的 最 大 距 


AWET ° 


这 样 ， 不 论 最 长 时 间 还 是 最 短 时 间 ， 都 只 要 对 每 只 蚂蚁 检查 一 次 就 好 
了 ， 这 是 O(n ) 时 间 的 算法 。 


对 于 限制 条 件 n < 106， 这 个 算法 是 够 用 的 ， 于 是 问题 得 解 。 
/ 输入 


int L, n; 

int x[MAX_N]; 

void solve() { 

/ 计算 最 短 时 间 

int minT = 0; 

for (int i = 0; i < n; i++) { 

minT = max(minT, min(x[i], L - x[i])); 
} 

/ 计算 最 长 时 间 

int maxT = 0; 

for (int i = 0; i < n; i++) { 

maxT = max(maxT, max(x[i], L - x[i])); 
} 


printf("%d %d\n", minT, maxT); 


} 


这 个 问题 可 以 说 是 考察 想象 力 类 型 问题 的 经 典 例子 。 有 很 多 这 样 的 问 
题 ， 虽 然 开始 不 太 明 日 ， 但 


想 通 之 后 ， 最 后 的 程序 却 是 出 乎 意料 地 简单 。 
1.6.3 难度 增加 的 抽签 问题 


如 果 将 最 开始 的 抽签 问题 中 关于 n 的 限制 条 件 改 为 1 < n < 1000 (题目 
描述 参见 1.1T) ， 那 么 应 


该 如 何 求解 呢 ? 最 初 的 四 重 循 环 算 法 是 O(n 4) 时 间 的 ， 将 n =1000 带 
A n 4 得 到 1012， 这 是 远 远 不 够 的 ， 


必须 改进 算法 。 


for (int a = 0; a < n; a++) { 

for (int b = 0; b < n; b++) { 

for (int c = 0; c < n; c++) { 

for (int d = 0; d < n; d++) { 

if (k[a] + k[b] + k[c] + k[d] == m) { 
f = true; 

} 

} 

} 
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上 面 是 最 初 所 记载 的 程序 的 循环 部 分 。 最 内 侧 关 于 d 的 循环 所 做 的 事 


就 是 


1 


检查 是 否 有 d 使 得 kDODkDkDkDma 
b 


C 
d 
通过 对 式 子 移 同 ， 束 能 得 到 为 一 种 表达 方式 


检查 是 否 有 d 使 得 kDmDkODkDkd 
a 


b 


C 
Mei, RAA k 中 所 有 元 妹 ， 判 断 是 否 有 mUkoDUkbDkc 。 
3 


证 我 们 着 眼 这 一 点 来 考虑 快速 的 检查 方法 。 虽 然 也 有 利用 数据 结构 优 
化 的 方法 (在 2.4 市 介绍 ) ， 


这 里 要 介绍 的 是 名 为 二 分 搜索 的 算法 。 
1. 二 分 搜索 与 O(n 3log n ) 的 算法 
4 


记 所 要 查找 的 值 mkaDKkbDkc 为 x。 预 先 把 数组 k 排 好 序 ， 然 后 看 
PRA, HA 


0D 如 果 它 比 x 小 ，x 只 可 能 在 
口 如 果 它 比 x 大 ，x 只 可 能 在 
5 


如 条 再 将 上 述 方 法 运用 在 已 经 减 半 的 x 的 存在 区 间 上 ，x 的 存在 区 间 
WLIO T ARR VA HER 


操作 束 可 以 不 断 缩 小 x ETEK IA), ERA UBER FRET > 


6 


7 
8 
9 
从 数列 中 查找 55 的 例子 


二 分 搜索 算法 每 次 将 候选 区 间 减 小 至 大 约 原来 的 一 半 。 因 此 ， 要 判断 
KENK n 的 有 序数 组 k 中 是 否 


10 


包含 x ， 只 要 反复 执行 约 log2 n 次 就 完成 了 。 二 分 查找 的 复杂 度 是 O 
dogmn) 时 间 的 ， 我 们 称 这 种 阶 的 


行 时 间 为 对 数 时 间 。 即 便 n 变 得 很 大 时 ， 对 数 时 间 的 算法 依然 非常 


速 。 


Ja 
天 


d 


11 


D MRTA TRACE IER, EA PRAF, ARM TA AOL 
的 某 边 的 值 。 
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即使 n 变 得 很 大 ，log2 n 也 依然 很 小 


n 


1 10 100 

1000 

106 109 

log2 n 0 

3 

7 

10 20 30 

将 最 内 侧 的 循环 玲 换 成 二 分 搜索 算法 之 后 ， 束 变 成 


O REF O (n log nm) 时间 
DER O (n Blog n ) 时 间 


n 3logn 比 nlogn 大 ， 所 以 这 里 合 起 来 当 作 O (n 3log n ) 时 间 
我 们 得 到 了 在 O(n 3log n ) 时 间 内 解决 


抽签 问题 的 算法 。 

/输入 

int n, m, KIMAX_N]; 

bool binary_search(int x) { 

1/ x 的 存在 范围 是 k[1], k[1+1), ..., k[r-1]. 
intl=0,r=n; 

/反复 操作 直到 存在 范围 为 至 
while (r -1 >= 1) { 
inti=(1+1/2; 

if (k[i] == x) return true; // 找到 x 
else if (k[i] <x)1=i + 1; 

else r = i; 

} 

/ 没 找到 x 

return false; 

} 


void solve() { 


°F 


H 
E, 


// 为 了 执行 二 分 查找 需要 先 排序 
sort(k, k + n); 

bool f = false; 

for (int a = 0; a < n; a++) { 

for (int b = 0; b < n; b++) { 

for (int c = 0; c < n; c++) { 

/将 最 内 侧 的 循环 替换 成 二 分 查找 
if (binary_search(m - k[a] - k[b] - k[c])) { 
f = true; 

} 

} 

} 

} 

if (£) puts("Yes"); 

else puts("No"); 

} 
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事实 上 ， 像 binary_search 这 样 的 函数 ， 多 数 情况 下 无 需 目 己 实 现 ， 可 以 
使 用 标准 实现 。 例 如 


1 
C++ 的 STL 中 就 含有 提供 基本 同样 功能 的 范 数 。 


2. O(n 2log n ) 的 算法 
2 


但 是 ， 将 n=1000 市 入 n 3log n ， 会 发 现 这 依然 是 远 远 不 够 的 ， 必 须要 
对 算法 做 进一步 优化 。 刚 才 我 


们 只 着 眼 于 四 重 循环 程序 中 最 内 层 的 循环 。 接 下 来 ， 让 我 们 着 上 腿 于 内 
侧 的 两 个 循环 。 


同 刚 才 一 样 的 思路 ， 内 侧 的 两 个 循环 是 在 


检查 是 否 有 c 和 d 侠 得 KDkDmDkDk。 
C 
d 
d 
b 


这 种 情况 并 不 能 直接 使 用 二 分 搜索 。 但 是 ， 如 果 预 完 枚 举 出 kc+ kd 所 
得 的 n2 个 数字 并 排 好 序 ， 便 可 


以 利用 二 分 搜索 了 GD。 
4 

该 算法 

O 排序 (n 2log n ) 时 间 
5 


O 循环 oO(n2logn) 时 间 


忆 的 也 十 O(n 2log n ) 时 间 。 这 样 束 可 以 确信 即便 n=1000 也 


对 了 。 

/ 输入 

6 

int n, m, k[MAX_N]; 

/ 保存 2 个 数 的 和 的 数列 
int kk[MAX_N * MAX_N]; 
7 


bool binary_search(int x) { 


// x 的 存在 范围 是 kk[1], kk[1+1], ..., kk[r-1]. 


intl=0,r=n*n; 
// 反复 操作 直到 存在 范围 为 空 


8 


while (r -1 >= 1) { 
inti=(1+r)/2; 

if (kk[i] == x) return true; // 找到 x 
else if (kk[i] < x) 1 =i + 1; 

else r = i; 

} 

9 

/ 没 找到 x 


台 巴 V7 


HE 


EM 


return false; 

} 

10 

void solveO { 

1/ 枚 举 k[c]+k[d] 的 和 


for (int c = 0; c < n; c++) { 


11 


O 确切 地 说 ， 去 除 重 复 后 n(n +1)/2 个 数字 就 够 了 ， 不 过 为 了 方便 可 以 
写成 榴 举 n 2 个 数字 。 
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for (int d = 0; d < n; d++) { 


kk[c * n + d] = k[c] + k[d]; 

} 

} 

/排序 以 便 进 行 二 分 搜索 

sort(kk, kk + n * n); 

bool f = false; 

for (int a = 0; a < n; a++) { 

for (int b = 0; b < n; b++) { 

1/ 将 内 侧 的 两 个 循环 蕉 换 成 二 分 搜索 


if (binary_search(m - k[a] - k[b])) { 
f = true; 

} 

} 

} 

if (£) puts("Yes"); 

else puts("No"); 

} 


本 问题 既 需要 二 分 搜索 这 一 基础 算法 知识 ， 也 需要 将 四 个 数 分 成 两 两 
考虑 的 想象 力 。 此 外 ， 像 这 


样 从 复杂 度 较 高 的 算法 出 发 ， 不 断 降 低 复杂 度 直到 满足 问题 要 求 的 过 
程 ， 也 是 设计 算法 时 常会 经 


历 的 一 个 过 程 。 
2 a 


DAA Dei 
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2.1 BAHN FRE” 


穷 竟 搜索 是 将 所 有 的 可 能 性 罗列 出 来 ， 在 其 中 寻找 答案 的 方法 。 这 里 
我 们 主要 介绍 深度 优先 搜 


索 和 广度 优先 搜索 这 两 种 方法 。 
2.1.1 递归 函数 


在 一 个 男 数 中 再 次 调用 该 函数 目 喘 的 行为 叫做 递归 ， 这 样 的 范 数 被 称 
(EVER o BE, BATA 


要 编写 一 个 计算 阶乘 的 函数 int fact(int n)， 当 然 ， 用 循环 来 实现 也 是 可 
以 的 。 但 是 根据 阶 


乘 的 递 推 式 n1= n x(n 一 1)!， 我 们 可 以 写成 如 下 形式 : 


int fact(int n) { 


if (n == 0) return 1; 
return n * fact(n - 1); 
} 


ERS NER, ER ACA RIERA ED EY ° CERI 
子 中 ， 当 n =0 时 fact 并 不 是 


继续 调用 目 身 ， 而 是 直接 返回 1。 如 采 没 有 这 一 条 件 存在 ， 函 数 束 会 无 
限 地 递归 下 去 ， 程 序 就 会 


RART ° 
fact 递 归 的 过 程 


我 们 再 来 试 试 编写 计算 斐 波 那 契 数列 的 函数 int fib(int n) ° RARR 
列 的 定义 是 a0=0、a1=1 


以 及 an= an [1+ anD02(n>1D。 这 里 ， 初 项 的 条 件 承 对 应 了 递归 的 终止 
条 件 。 数 列 的 定义 直接 写成 贸 数 


RAL LAT ° 


2.1 最 基础 的 “ 穷 网 搜索 ”27 


int fib(int n) { 


1 

if (n <= 1) return n; 

return fib(n - 1) + fib(n - 2); 
} 


实际 使 用 这 个 函数 时 ， 即 使 是 求 fib(40) 这 样 的 n BUNTAR, BEE 
费 相当 长 的 时 间 。 这 是 


2 
因为 这 个 函数 在 递归 时 ， 会 像 下 图 一 样 按照 指数 级 别 扩展 开 来 。 
3 


4 


fib(10) 递归 的 过 程 


在 斐 波 那 契 数列 中 ， 如 果 fibAo) 的 是 一 定 的 ， 无 论 多 少 次 调用 都 会 得 
到 同样 的 结 末 。 因 此 如 采 


5 


计算 一 次 之 后 ， 用 数列 将 结果 存储 起 来 ， 便 可 优化 之 后 的 计算 。 (上 
图 中 ) fib(10) 被 调用 时 同 


样 的 n 被 计算 了 很 多 次 ， 因 此 可 以 获得 很 大 的 优化 空间 。 这 种 方法 是 
出 于 记忆 化 搜索 或 者 动态 规 


划 的 想法 ， 之 后 我 们 会 介绍 。 
6 

int memo[MAX_N + 1]; 

int fib(int n) { 

if (n <= 1) return n; 

7 

if (memo[n] != 0) return memo[n]; 
return memo[n] = fib(n - 1) + fib(n - 2); 
} 

8 

2.1.2 栈 


栈 (Stack) 是 支持 push 和 pop 两 种 操作 的 数据 结构 。push 是 在 栈 的 顶端 
放 入 一 组 数据 的 操作 。 反 


之 ，pop 走 从 其 顶端 取出 一 组 数据 的 操作 。 因 此 ， 最 后 进入 栈 的 一 组 效 
据 可 以 最 先 被 取出 (这 种 


9 
行为 被 叫做 LIFO: Last In First Out， 即 后 进 先 出 ) 。 


通过 使 用 数组 或 者 列表 等 结构 可 以 很 容易 实现 栈 ， 不 过 C++、Java 等 程 
序 语 言 的 标准 库 已 经 为 我 


们 准备 好 了 这 一 常用 结构 ， 在 比赛 中 需要 时 不 妨 使 用 它们 。C++ 的 标准 
库 中 ，stack::pop 完 成 


10 


的 仅仅 是 移 除 最 顶端 的 数据 。 如 果 要 访问 最 顶端 的 数据 ， 需 要 使 用 
stack::top 函 数 (这 个 操作 


通常 也 被 称 为 peek) 。 


11 


| i 
vs) || pi — 
JE 
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=| — 
ie 
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六 数 调 用 的 过 程 是 通过 使 用 栈 实现 的 。 因 此 ， 尼 归 函 数 的 递归 过 程 也 
可 以 改 用 栈 上 的 操作 来 实现 。 


现实 中 需要 如 此 改写 的 场合 并 不 多 ， 不 过 作为 使 用 栈 的 练习 试 试看 也 
征 不 铺 的 。 以 下 是 使 用 


stack 的 例子 : 

#include <stack> 

#include <cstdio> 

using namespace std; 

int main() { 

stack<int> s; // 声明 存储 int 类 型 数据 的 栈 
s.push(1); / {} > {1} 

s.push(2); // {1} > {1,2} 

s.push(3); // {1,2} > {1,2,3} 
printf("%d\n", s.top(); // 3 

s.pop(); // 从 栈 顶 移 除 {1,2,3} > {1,2} 
printf("%d\n", s.top()); // 2 

s.pop(); // {1,2} > {1} 

printf("%d\n", s.top()); // 1 

s.pop(); // {1} > € 

return 0; 


} 


2.1.3 队列 


队列 (Queue) 与 栈 一 样 文 持 push 和 pop 两 个 操作 。 但 与 栈 不 同 的 是 ， 
pop 完 成 的 不 是 取出 最 顶端 


的 元 素 ， 而 是 取出 最 瓜 剖 的 元 素 。 也 束 古 说 最 初 放 入 的 元 素 能 够 最 先 


被 取出 (这 种 行为 被 叫做 FIFO: 
First In First Out， 即 先进 先 出 ) 。 


¡MI u 
vs) | |p — 


E 


My (MS 
an 
ll 


y TT A 


3 
4 
5 
队列 的 操作 


如 同 栈 一 样 ，C++、Java 等 的 标准 库 也 预 置 了 队列 。Java 与 C++ 中 的 男 
数 的 名 称 与 用 途 稍 有 不 同 ， 


因此 使 用 时 要 注意 。 此 外 ， 在 C++ 中 queue::front 是 用 来 访问 最 底 端 数据 
的 函数 。 以 下 是 使 用 


6 

queue 的 例子 : 

#include <queue> 

#include <cstdio> 

7 

using namespace std; 

int main() { 

queue<int> que; // 声明 存储 int 类 型 数据 的 队列 
que.push(1); // {} > {1} 

8 

que.push(2); // {1} > {1,2} 


que.push(3); // {1,2} > {1,2,3} 


printf("%d\n", que.frontO); // 1 
que.pop(); // 从 队 尾 移 除 {1,2,3} > {2,3} 
printf("%d\n", que.front()); // 2 
9 

que.pop(); // {2,3} > {3} 
printf("%d\n", que.front()); // 3 
que.pop(); // {3} > {} 

return 0; 

} 

10 

2.1.4 深度 优先 搜索 

11 


深度 优先 搜索 (DFS, Depth-First Search) 是 搜索 的 手段 之 一 。 它 从 某 
个 状态 开始 ， 不 断 地 转移 
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状态 直到 无 法 转移 ， 然 后 回 退 到 前 一 步 的 状态 ， 继 续 转 移 到 其 他 状 
仿 ， 如 此 不 断 重复 ， 直 至 找到 


最 终 的 解 。 例 如 求解 数 独 ， 站 先 在 菏 个 格子 内 填 入 适当 的 数 子 ， 然 后 
再 继续 在 下 一 个 格子 内 填 入 


数字 ， 如 此 继续 下 去 。 如 采 发 现 某 个 格子 无 解 了 ， 吏 放弃 前 一 个 格子 
上 选择 的 数字 ， 改 用 其 他 可 


行 的 数字 。 根 据 深 度 优 先 搜索 的 特点 ， 采 用 递归 函数 实现 比较 简单 。 
状态 转移 的 顺序 


我 们 来 试 痢 解答 一 下 下 面 的 题目 : 
部 分 和 问题 


给 定 整 数 a1、a2、...、an ， 判 断 是 否 可 以 从 中 选 出 奋 干 数 ， 使 它们 
的 和 恰好 为 k。 


限制 条 件 
[1<nz<20 

O 0108 < ai < 108 
O [108 < k < 108 
样 例 1 

输入 

n=4 

a={1,2,4,7} 

k=13 

输出 


Yes (13=2+4+7) 


2.1 最 基础 的 “ 穷 奖 搜索 ”31 
样 例 2 

1 

输入 

n=4 
a={1,2,4,7} 
k=15 

2 

输出 

No 

3 


从 a 1 开始 按 顺 序 决定 每 个 数 加 或 不 加 ， 在 全 部 mn 个 数 都 决定 后 再 判断 
它们 的 和 是 不 是 k 即 可 。 因 为 


Se 


状态 数 是 2 n +1， 所 以 复杂 度 是 O (2 n)。 如 何 实现 这 个 搜索 ， 请 参见 
下 面 的 代码 。 注 意 a 的 下 标 与 题 


目 描述 中 的 下 标 偶 移 了 1。 在 程序 中 使 用 的 是 0 起 始 的 下 标 规 则 ， 题 目 
描述 中 则 是 1 开始 的 ， 这 一 


4 
点 要 注意 避免 搞 混 。 


5 


6 

7 

8 

状态 转移 的 样子 

/输入 

int a[MAX_N]; 

9 

int n, k; 

/ 已 经 从 前 ji 项 得 到 了 和 sum， 然 后 对 于 i 项 之 后 的 进行 分 文 
bool dfs(int i, int sum) { 

/ 如果 前 n 项 都 计算 过 了 ， 则 返回 sum 是 否 与 k 相 等 
10 

if (i == n) return sum == 


/不 加 上 afi] 的 情况 


if (dfs(i + 1, sum)) return true; 
11 
/加 上 ai 的 情况 
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if (dfs(i + 1, sum + ali])) return true; 

// 无 论 是 否 加 上 a 都 不 能 竣 成 k 环 返回 false 
return false; 

} 

void solve() { 

if (dfs(0, 0)) printf("Yes\n"); 

else printf("No\n"); 

} 


深度 优先 搜索 从 最 开始 的 状态 出 发 ， 避 历 所 有 可 以 到 达 的 状态 。 由 此 
可 以 对 所 有 的 状态 进行 操作 ， 


或 者 列举 出 所 有 的 状态 。 


Lake Counting (POJ No.2386) 


有 一 个 大 小 为 N x M 的 园子 ， 雨 后 积 起 了 水 。 八 连通 的 积 水 补 认为 十 
连接 在 一 起 的 。 请 求 出 


园子 里 总 共有 多 少 水 洼 ? 《 八 连通 指 的 是 下 图 中 相对 W 的 * 的 部 分 ) 


KK 


*W* 


限制 条 件 

ON,M<100 

样 例 

输入 

N=10, M=12 

园子 如 下 图 (WEB, ERAK) 


.WAN 
输出 

3 

2.1 最 基础 的 “ 穷 奖 搜索 ”33 


从 任意 的 W 开 始 ， 不 集 地 把 邻接 的 部 分 用 "代替 。1 次 DFS 后 与 初始 的 


这 个 W 连 接 的 所 有 W 束 部 被 叔 


1 


换 成 了 '"， 因 此 直到 图 中 不 再 存在 W 为 上 上， 


案 了 。8 个 方 同 共 对 应 了 8 种 


状态 转移 ， 每 个 格子 作为 DFS 的 参数 至 多 被 调用 一 次 ， 所 以 复杂 度 为 


(8xNxM)=O(NxM)。 
/ 输入 

2 

int N, M; 

char field[MAX_N][MAX_M + 1]; // 园子 
/ 现在 位 置 (x,y) 

void dfs(int x, int y) { 

3 

I 将 现在 所 在 位 置 奉 换 为 . 
field[x]Ly] ="; 
NENA 


for (int dx = -1; dx <= 1; dx++) { 


总 共 进 行 DFS 的 次 数 就 是 答 


O 


for (int dy = -1; dy <= 1; dy++) { 

4 

// 问 X 方 向 移动 dk， 和 疝 y 方 向 移动 dy， 移 动 的 结果 为 (nx,ny) 
int nx = x + dx, ny = y + dy; 

/ 判断 (nx,ny) 是 不 是 在 园子 内 ， 以 及 是 否 有 积 水 


if (0 <= nx && nx < N && 0 <= ny && ny < M && field[nx][ny] == 'W') 
dfs(nx, ny); 


} 
5 


} 


return ; 

} 

void solve() { 

6 

int res = 0; 

for (int i = 0; i < N; i++) { 
for (int j = 0; j < M; j++) { 
if (field[i][j] == 'W') { 

/ 从 有 W 的 地 方 开始 dfs 
7 


dfs(i, j); 


res++; 


printf("%d\n", res); 


8 
} 
2.1.5 宽度 优先 搜索 
9 


宽度 优先 搜索 (BES, Breadth-First Search) 也 是 搜索 的 手段 之 一 。 它 
与 深度 优先 搜索 类 似 ， 从 某 


个 状态 出 发 探索 所 有 可 以 到 达 的 状态 。 


与 深度 优先 搜索 的 不 同 之 处 在 于 搜索 的 顺序 ， 宽 度 优先 搜索 总 是 先 搜 
索 距离 初始 状态 近 的 状态 。 


10 


也 就 是 说 ， 它 是 按照 开始 状态 ~ 只 需 1 次 转移 就 可 以 到 达 的 所 有 状态 ~ 
AREA 以 到 达 的 


所 有 状态 二 oo. 这 样 的 顺序 进行 搜索 。 对 于 同一 个 状态 ， 宽 度 优先 搜 
索 只 经 过 一 次 ， 因 此 复杂 度 


为 O (状态 数 x 转 移 的 方式 )。 


11 


V V 
De 
RUN 
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状态 转移 的 顺序 


深度 优先 搜索 ( 隐 式 地 ) 利用 了 栈 进行 计算 ， 而 宽度 优先 搜索 则 利用 
了 队列 。 搜 索 时 首先 将 初始 


状态 添加 到 队列 里 ， 此 后 从 队列 的 最 前 闻 不 断 取 出 状态 ， 把 从 该 状态 
可 以 转移 到 的 状态 中 尚未 访 


问 过 的 部 分 加 入 队列 ， 如 此 往复 ， 直 至 队列 被 取 空 或 找到 了 问题 的 
解 。 通 过 观察 这 个 队列 ， 我 们 


00 a 


TAS A Ce 


给 定 一 个 大 小 为 WNx MINE E EIR A er BE, Be YD 
向 邻接 的 上 下 左右 四 格 


的 通道 移动 。 请 求 出 从 起 点 到 终点 所 需 的 最 小 步 数 。 请 注意 ， 本 题 假 
定 从 起 点 一 定 可 以 移动 


到 终点 。 
限制 条 件 
ON,M<100 
样 例 

输入 


N=10, M=10 (迷宫 如 下 图 所 示 。' 太 ，'"，'S'，'G' 分 别 表示 墙壁 、 通 道 、 


HSHHHHHHH 


wens #..H 


HAH HH 


Mri 


HHHH HHHH 


nH.. H 


HEHEHE H 


vc ios 


HHRH HHH. 


GA 


2.1 最 基础 的 “ 穷 奖 搜索 ”35 
输出 

1 

22 


宽度 优先 搜索 按照 距 开始 状态 由 近 及 远 的 顺序 进行 搜索 ， 因 此 可 以 很 
容易 地 用 来 求 最 短路 径 、 最 


少 操 作 之 类 问题 的 答案 。 这 个 问题 中 ， 状 态 仅 仅 古 目前 所 在 位 置 的 坐 
标 ， 因 此 可 以 构造 成 pair 


2 


或 者 编码 成 int 来 表达 状态 。 当 状态 更 加 复杂 时 ， 就 需要 封装 成 一 个 类 
来 表示 状态 了 。 转 移 的 方 


式 为 四 方向 移动 ， 状 态 数 与 迷宫 的 大 小 是 相等 的 ， 所 以 复杂 度 是 O (4x 
NxM)=O(NXM)° 


宽度 优先 搜索 中 ， 只 要 将 已 经 访问 过 的 状态 用 标记 管理 起 来 ， 就 可 以 
很 好 地 做 到 由 近 及 远 的 搜索 。3 


这 个 问题 中 由 于 要 求 最 短 距离 ， 不 妨 用 d[N][M] 数 组 把 最 短 距离 保存 起 
来 。 初 始 时 用 充分 大 的 第 


数 INF 来 初始 化 它 ， 这 样 疝 未 到 达 的 位 置 束 古 INF， 也 束 同 时 起 到 了 标 
记 的 作用 。 


4 


里 然 到 达 终 点 时 就 会 停止 搜索 ， 可 如 采 继 续 下 去 直到 队列 为 空 的 话 ， 
束 可 以 计算 出 到 各 个 位 置 的 


最 短 距离 。 此 外 ， 如 果 搜 索 到 最 后 ，d 依 然 为 INF 的 话 ， 便 可 得 知 这 个 
位 置 束 是 无 法 从 起 点 到 达 的 


WE > 
5 


在 今后 的 程序 中 ， 使 用 像 INF 这 样 充 分 大 的 常数 的 情况 还 很 多 。 不 把 
INE 当 作 例 外 ， 而 是 直接 参 


普通 运算 的 情况 也 很 常见 。 这 种 情况 下 ， 如 果 INEF 过 大 就 可 能 带 来 溢 
出 的 危险 。 


假设 INF=231-1。 例 如 想 用 d[nx][ny]=min(d[nx][ny], d[x][y]+1D) 来 更 新 
d[nx][ny], MZ 


6 


发 生 INF+1=-231 的 情况 。 这 一 问题 中 d[x][y] 总 不 等 于 INE， 所 以 没有 问 
题 。 但 是 为 了 防止 这 样 


的 问题 ， 一 般 会 将 INEF 设 为 放大 2~4 倍 也 不 会 溢出 的 大 小 〈 可 参考 2.5 人 T 
Floyd-Warshall 算 法 等 ) 。 


因为 要 癌 4 个 不 同方 同 移动 ， 用 dx[4] 和 dy[4 两 个 数组 来 表示 四 个 方 同 问 


量 。 这 样 通过 一 个 循 

7 

环 殉 可 以 实现 四 方 癌 移动 的 忆 历 。 

const int INF = 100000000; 

/ 使 用 pair 表 示 状 态 时 ， 使 用 typedef 会 更 加 方便 一 些 
8 

typedef pair<int, int> P; 

/输入 


char maze[MAX_N][MAX_M + 1]; V 表示 迷宫 的 字符 串 的 数组 


int N, M; 

9 

int sx, sy; // 起 点 坐标 

int gx, gy; // 终点 坐标 

int dIMAX_N]IMAX_M]; // 到 各 个 位 置 的 最 短 距离 的 数组 
10 

1 ANT HSS E 

int dx[4] = {1, 0, -1, 0}, dy[4] = {0, 1, 0, -1}; 

/ 求 从 (sx, sy) 到 (gx, gy) 的 最 短 距离 

/ 如 果 无 法 到 达 ， 则 是 INF 

int bfs() { 

11 
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queue<P> que; 

/把 所 有 的 位 置 都 初始 化 为 INF 

for (int i = 0; i < N; i++) 

for (int j = 0; j < M; j++) d[illj] = INF; 

/ 将 起 点 加 入 队列 ， 并 把 这 一 地 点 的 距离 设置 为 0 


que.push(P(sx, sy)); 


d[sx][sy] = 0; 


/ 不 断 循环 直到 队列 的 长 度 为 0 

while (que.size()) { 

/ 从 队列 的 最 前 端 取出 元 素 

P p = que.front(); que.pop(); 

/如 果 取 出 的 状态 已 经 是 终点 ， 则 结束 搜索 
if (p.first == gx && p.second == gy) break; 

1/ 四 个 方 同 的 循环 

for (int i = 0; i < 4; i++) { 

// 移动 之 后 的 位 置 记 为 (nx, ny) 

int nx = p.first + dx[i], ny = p.second + dy[i]; 


/判断 是 否 可 以 移动 以 及 是 否 已 经 访问 过 (d[nx]j[ny]!=INF 即 为 已 经 访 
问 过 ) 

if (0 <= nx && nx < N && 0 <= ny && ny < M && maze[nx][ny] != '#' 
&& d[nx][ny] == INF) { 


/ 可 以 移动 的 话 ， 则 加 入 到 队列 ， 并 且 到 该 位 置 的 距离 确定 为 到 p 的 距 
ray +1 


que.push(P(nx, ny)); 

d[nx][ny] = d[p.first][p.second] + 1; 
} 

} 

} 


return d[gx][gy]; 


} 

void solve() { 

int res = bfs(); 
printf("%d\n", res); 
} 


ee BE, SEE REE EAI 
， 因 此 需要 对 所 有 状态 进行 


un, 度 优 先 搜索 也 是 可 以 的 。 但 是 递归 函数 可 以 很 简短 地 编 
， 而 且 状 态 的 管理 也 更 们 


单 ， 所 以 大 多 数 情况 下 还 是 用 深度 优先 搜索 实现 。 反 之 ， 在 求 取 最 短 
路 时 深度 优先 搜索 需要 反复 


过 同样 的 状态 ， 所 以 此 时 还 是 使 用 宽度 优先 搜索 为 好 。 


宽度 优先 搜索 会 把 状态 逐个 加 入 队列 ， 因 此 通常 需要 与 状态 数 成 正比 
的 内 存 空 ao RZ, RER 


先 搜 索 是 与 最 大 的 递归 深度 成 正比 的 。 一 般 与 状态 数 相 比 ， 递 归 的 深 
度 并 不 会 太 大 ， 所 以 可 以 认 


为 深度 优先 搜索 更 加 节省 内 存 。 


此 外 ， 也 有 采用 与 宽度 优先 搜索 类 似 的 状态 转移 顺序 ， 并 且 注 重 节 约 
内 存 占 用 的 迭代 加 深 深度 优 


先 搜索 (IDDFS, Iterative Deepening Depth-First Search) 。IDDFS 是 一 
种 在 最 开始 将 深度 优先 搜索 


的 递归 次 数 限 制 在 1 次 ， 在 找到 解 之 前 不 断 增 加 递归 深度 的 方法 。 这 种 
方法 会 在 4.5 市 详细 介绍 。 


2.1 最 基础 的 “ 穷 济 搜索 ”37 


1 
2.1.6 特殊 状态 的 枚 举 


里 然 生成 可 行 解 空间 多 数 采用 深度 优先 搜索 ， 但 在 状态 空间 比较 特殊 
时 其 实 可 以 很 商 短 地 实现 。 


比如 ，C++ 的 标准 库 中 提供 了 next_permutation 这 一 函数 ， 可 以 把 了 个 元 
RE n ! 种 不 同 的 排列 生 


2 


o 又 或 者 ， 通 过 使 用 位 运算 ， 可 以 枚 举 从 nn 个 元 于 中 取出 kk 个 
JÆ k 


C 种 状态 或 是 菜 个 集合 


n 
中 的 全 部 子 集 等 。3.2 节 将 介绍 如 何 利用 位 运算 枚 举 状态 。 
bool used[MAX_N]; 

int perm[MAX_N]; 

3 

/ 生成 {0,1,2,3,4,…,n-1} 的 n! 种 排列 

void permutation1(int pos, int n) { 

if (pos == n) { 

4 

po 

* 这 里 编写 需要 对 perm 进 行 的 操作 


el 


return ; 
} 

5 

/ 针对 perm 的 第 pos 个 位 置 ， 究 竟 使 用 0~n-1 中 的 哪 一 个 进行 循环 
for (int i = 0; i < n; i++) { 

if (lused[i]) { 

perm[pos] = i; 

// 已 经 被 使 用 了 ， 所 以 把 标志 设置 为 true 

6 

used[i] = true; 

permutation1(pos + 1, n); 

/返回 之 后 把 标志 复位 

used[i] = false; 

} 


} 
7 


return ; 
} 
#include <algorithm> 


8 


/即使 有 重复 的 元 素 也 会 生成 所 有 的 排列 
// next_permutation 是 按照 字典 序 来 生成 下 一 个 排列 的 


int perm2[MA X_N]; 


void permutation2(int n) { 

9 

for (int i = 0; i < n; i++) { 

perm2[i] = i; 

} 

do { 

Pr 

*+ 这 里 编写 需要 对 perm2 进 行 的 操作 

10 

El 

} while (next_permutation(perm2, perm2 + n)); 
// 所 有 的 排列 都 生成 后 ，next_permutation 会 返回 false 
return ; 

} 

11 


V V 
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ye 


38 第 2 音 初 出 茅 庐 一 一 初级 篇 
2.1.7 BAX 


MEEL, FRAT ATA AREA, AAA 
时 ， 复 杂 度 也 会 相应 变 大 。 


比如 个 元 素 进 行 排列 时 状态 数 总 共有 个， 复杂 度 也 就 成 了 O(n 
!)。 这 样 的 话 ， 即 使 n=15 计 算 也 


很 难 较 早 终止 。 这 里 简单 介绍 一 下 此 类 情形 要 如 何 进行 优化 。 


深度 优先 搜索 时 ， 有 时 早已 很 明确 地 知道 从 当前 状态 无 论 如 何 转 移 都 
` 会 存在 解 。 这 种 情况 下 ， 


不 再 继续 搜索 而 是 直接 跳 过 ， 这 一 方法 个 称 作 静 枝 。 


我 们 回想 一 下 深度 优先 搜索 的 例题 “部 分 和 问题 ”。 这 个 问题 中 的 限制 
条 件 如 果 变 为 0< ai <108, 


那么 在 递归 中 只 要 sum 超 过 k 了 ， 此 后 无 论 选 择 哪些 数 都 不 可 能 让 sum 等 
于 k， 所 以 此 后 没有 必要 


继续 搜索 。 


BIRIT OL 
关于 更 多 更 高 级 的 搜索 手段 ， 我 们 会 在 4.5 节 进行 详细 介绍 
专栏 栈 内 存 和 堆 内 存 


调用 函数 时 ， 主 调 的 函数 所 拥有 的 局 部 变量 等 信息 需要 存储 在 特定 的 
内 存 区 域 。 这 个 区 域 被 


称 作 栈 内 存 区 。 另 一 方面 ， 利 用 new 或 者 malloc 进行 分 配 的 内 存 区 域 
被 称 为 堆 内 存 。 


栈 内 存在 程序 启动 时 被 统一 分 配 ， 此 后 不 能 再 扩大 。 由 于 这 一 区 域 有 
上 限 ， 所 以 函数 的 递归 


深度 也 有 上 限 。 虽 然 与 函数 中 定义 的 局 部 变量 的 数目 有 关 ， 不 过 一 般 
情况 下 C 和 C++ 中 进行 


上 万 次 的 递归 是 可 以 的 。 在 Java 中 ， 在 执行 程序 时 可 以 用 参数 指定 栈 
的 大 小 。 不 同 的 程序 设 


计 竞赛 所 采用 的 设置 各 有 不 同 ， 建 议 大 家 预先 进行 确认 。GCJ 的 话 ， 
程序 是 在 自己 的 机 器 上 


执行 的 ， 所 以 可 以 自行 设置 参数 。 


全 局 变量 被 保存 在 扒 内 存 区 。 通 常 不 推荐 使 用 全 局 变量 ， 但 是 在 程序 
RUTA, HT REN 
通 单 不 是 那么 多 ， 并 且 香 负 会 有 多 个 函数 访问 同一 个 数组 ， 因 此 利用 
Sipe BRA te > HE 
bs ENANA ECAC, IBERIA CIEE, RG 
在 堆 内 存 上 可 以 减少 栈 


洲 出 的 危险 。 同 时 ， 通 第 只 需 定义 满足 最 大 需要 的 数列 大 小 ， 但 如 采 
N E 


Penn en 字符 串 末 尾 的 \0' 的 空间 之 类 的 漏 
洞 。 
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1 
2.7 一 起 来 挑战 GCJ 的 题目 (1) 

2 

让 我 们 一 起 运用 迄今 所 介绍 的 技巧 ， 实 际 挑战 一 下 GCJ 的 题目 吧 。 
3 


2.7.1 Minimum Scalar Product 
4 
Minimum Scalar Product (2008 Round1A A) 


AMANTE v 1=(x 1, x 2, ..., xn Al v 2=(y 1,y2,...,yn)， 人 允许 任意 交换 
v1 和 v2 各 目的 分 量 的 顺序 。 请 


计算 v1 和 v2 的 内 积 x1y 1+...+ xnyn 的 最 小 值 。 
5 
限制 条 件 


Small 


[1<n<8 


O 01000 < xi, yi < 1000 


6 


Large 


[ 100 <n < 800 


O 0100000 < xi , yi < 100000 


7 
样 例 1 

输入 

8 

n=3 

v= (1, 3, -5) 
1 

y =(-2, 4, 1) 
2 

输出 


9 


令 V=(-5, 1, 3),V= (4, 1, -2) 束 可 以 得 到 最 小 值 v x v = -25 


1 


2 


1 
2 

样 例 2 

10 

输入 

n=5 

v =(1, 2, 3, 4, 5) 

1 

v= (10; 1,0, 1) 

2 

11 
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输出 

4-v = (1, 2, 3, 4, 5), v=(1, 1, 1, 0, 0) 就 可 以 得 到 最 小 值 6 
1 

2 


v1 和 v2 各 自 的 分 量 的 顺序 都 可 以 任意 交换 ， 因 此 可 以 先 把 v1 的 顺序 
固定 下 来 只 交换 v 2 的 顺序 。 为 了 


方便 分 析 先 将 v 1 按 升序 排 好 序 。 接 下 来 枚 举 v 2 的 分 量 所 有 的 排列 顺 
序 ， 一 共有 了 ! 种 排列 ， 还 需要 


对 每 种 排列 计算 内 积 ， 总 的 复杂 度 是 O(n !x n )。 这 个 算法 在 Small 的 
情况 ， 因 为 n< 8 所 以 没有 问题 ， 


但 在 Large 的 情况 时 就 远 远 不 够 了 。 


在 此 隐约 会 觉得 把 y 2 按 降序 或 者 升序 排序 的 话 ， 所 得 的 内 积 是 最 小 
的 。 事 实 上 ， 如 果 将 v 2 按 降 序 


排序 的 话 ， 所 给 的 两 个 样 例 都 能 够 得 到 最 小 值 。 这 个 设想 的 确 是 对 
的 ，v1 和 v2 的 内 积 在 将 v 1 按 升 


序 ， 将 v2 按 降序 排序 时 取得 最 小 值 。 下 面 我 们 来 证 明 这 一 设想 。 首 先 
考 虚 n=2 的 情况 。 


考 虚 vy 1=(x 1,x 2), v2=(y1,y 2), 假设 v 1 已 经 按 升序 排 好 序 了 ， 即 有 
x1<x2, 比较 x1xy1l1+x2xy2 和 x2xy1l 


+XxX1xy2 的 大 小 关系 。 
x1xy1+x2xy2[]x2xy1[]x1xy2=x1(y1[l]y2)+x2(y2l]y 
1)=(x1 0x2) (y10y 2) BT, y1>y2[x1xy1+x2xy2<x2 
xy1+x1xy2。 因 此 n=2 时 结论 成 立 。 


接 下 来 考虑 n 大 于 2 的 情况 。 如 果 v2 不 是 按 降序 排序 的 ， 那 么 存在 i <j 
使 得 yi<yj， 根 据 对 n =2 的 情 


况 的 分 析 可 以 知道 ， 交 换 yi 入 后 束 得 到 了 更 小 的 内 积 。 因 此 ， 当 将 
v 2 按 降序 排序 时 ， 所 得 的 内 积 


最 小 。 


数组 排序 的 复杂 度 为 O(nlog n )， 所 以 发 现 这 一 结论 后 ， 只 要 做 两 次 
排序 就 可 以 简单 高 效 地 计算 管 


案 了 。 那 么 ， 要 怎样 才能 去 想到 这 个 设想 呢 ? 
首 允 就 是 靠 直 觉 。 比 较 容 易 想 到 的 姑且 移 排序 试 试看 。 


第 二 个 就 古 样 例 。 原 本 的 问题 揪 述 中 并 没有 说 明 v1 和 v2 在 什么 情况 
下 取得 最 小 值 。 但 是 ， 因 为 样 


例 的 规模 比较 小 ， 所 以 可 以 手工 验算 。 因 而 可 以 找到 使 得 内 积 最 小 的 v 
1 和 和 v2。 通过 观察 对 应 的 v1 


和 v2， 想 要 得 到 刚才 的 结论 也 并 不 古 那 么 难 。 


第 三 个 就 是 像 证 明 中 的 那样 ， 从 像 n=2 这 样 小 规模 的 情况 出 发 ， 推 广 
到 一 般 的 情况 从 而 证 明 结 


论 。 在 这 个 问题 中 ， 很 容易 证 明 对 n =2 的 情况 结论 成 立 ， 接 着 同样 可 
以 证 明 对 一 般 的 情况 结论 


AB BOIL ° 


最 后 需要 注意 的 是 ， 即 使 推导 出 了 正确 的 算法 ， 如 果 程 序 的 实现 中 有 
漏洞 的 话 也 十 徒 亡 。 在 这 道 


题 中 ，v1 的 分 量 和 v 2 的 分 量 的 乘积 可 能 会 导致 32 位 整数 盗 出。 即使 
认为 找到 了 正解 ， 不 到 最 后 所 


交 结 采 正 确 的 那 一 刻 ， 都 不 能 挥 以 轻 心 呢 。 


u 出 由 
AO 


= 0000 
0000000 
goo 
AA 


2.7 一 起 来 挑战 GCJ 的 题目 (1) 127 


typedef long long II; 

1 

/输入 

int n; 

int vIIMAX_N], v2[MAX_N]; 
2 

void solve() { 


sort(v1, v1 + n); 


sort(v2, v2 + n); 

ll ans = 0; 

for (int i = 0; i < n; i++) ans += (11)v1[i] * v2[n - i - 1]; 3 
printf("%lld\n", ans); 

} 

2.7.2 Crazy Rows 

4 

Crazy Rows (2009 Round2 A) 

5 


给 定 一 个 由 0 和 1 组 成 的 矩阵 。 只 人 允许 交换 相 邻 的 两 行 (第 i 行 和 第 i 
+147) ， 要 把 矩阵 化 成 


下 三 角 和 矩阵 〈 主 对 角 线 上 方 的 元 素 都 是 0) ， 最 少 需 要 交换 几 次 ?输入 
的 矩阵 保证 总 能 化 成 下 


= FA FEB 。 


6 


7 

行 交 换 
变换 后 
限制 条 件 
Small 


[1<N<8 


8 

Large 

[]}1<N<40 

9 

样 例 1 

输入 

N=2 

10 

JEBE 

10 

11 

11 

128 第 2 章 初出 茅 访 一 一 初级 篇 
输出 

0( 输 入 已 经 是 下 三 角 和 矩阵 ) 
样 例 2 

输入 

N=3 

JEBE 


001 


100 
010 

输出 

2( 交 换 第 1 行 和 第 2 行 ， 再 交换 第 2 行 和 第 3 行 ， 得 到 下 三 角 和 矩阵 ) 
样 例 3 

输入 

N=4 

和 矩阵 


1110 


1100 
1100 
1000 
输出 
4 


最 先 想 到 的 是 尝试 所 有 NN ! 种 交换 方案 。 但 在 Large 中 ， 由 于 最 大 的 N 
=40， 这 当然 是 行 不 通 的 。 
暂且 移 考 虑 一 下 最 后 应 该 把 哪 一 行 交 换 到 第 1 行 。 最 后 的 第 1 行 应 该 具 
有 00...0 或 是 10...0 的 形式 。 
可 以 交换 到 第 1 行 的 行当 然 也 可 以 交换 到 第 2 及 之 后 的 行 ， 当 有 多 个 满 
足 条 件 的 行 时 ， 选 择 离 第 1 


行 近 的 行 对 应 的 最 终 费 用 要 小 。 大 家 肯定 都 已 注意 到 了 这 一 点 吧 。 有 
兴趣 的 读 着 不 妨 目 己 证 明 


—Fo 


确定 第 1 行 之 后 ， 束 没有 必要 再 移 动 它 了 ， 于 是 对 于 之 后 的 行 束 可 以 以 
同样 的 思路 处 理 。 


在 这 道 题 中 ， 每 行 的 0 和 1 的 位 置 并 不 重要 ， 只 要 知道 每 行 最 后 一 个 1 所 
ENEE T ° WR 


先 将 这 些 位 置 预 先 计 算 好 ， 那 么 束 能 降低 行 交 换 时 的 复杂 上 度 。 直 接 按 
和 矩阵 的 形式 处 理 的 复杂 度 是 


O(N 3)， 而 预先 计算 后 再 处 理 的 复杂 度 降 为 O(ND > 
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// 输入 


1 


int N; 


int M[MAX_NJ[MAX_N]; // 矩阵 


int a[MAX_N]; // ali ar Bi TER AI E 10-1 
2 

void solve() { 

int res = 0; 

for (int i = 0; i < N; i++) { 

ali] = -1 / WR BÁTAR, M1 


for (int j = 0; j < N; j++) { 


3 
if (ML[LiJD] == 1) ali] = j; 
} 
} 


for (int i = 0; i < N; i++) { 


int pos = -1; // 要 移动 到 第 i 行 的 行 
for (int j = i; j < N; j++) { 
4 

if (alj] <= i) { 

pos = j; 

break; 

} 

i 

5 

/ 完成 交换 

for (int j = pos; j > i; j--) { 
swap(alj], alj - 1]); 
TeS 十 十 ; 

6 


} 


} 

printf("%d\n", res); 

} 

7 

2.7.3 Bribe the Prisoners 

8 

Bribe the Prisoners (2009 Round 1C C) 


如 下 图 所 示 ， 一 个 监狱 里 有 P 个 并 排 着 的 牢房 。 从 左 至 右 依次 编号 为 
1,2, ..., P。 最 初 所 有 的 牢 


房 里 都 住 着 一 个 囚犯。 相 邻 的 两 个 牢房 之 间 有 一 个 窗户 ， 可 以 通过 它 
与 相 邻 牢房 里 的 办 犯 对 话 。 


9 


10 
监狱 的 情况 


INE PEN ENA > MARC TES TAIZ, 其 相 邻 的 牢房 里 
的 内 犯 束 会 知道 ， 因 而 


11 


发 怒 骏 动 。 所 以 ， 释 放 某 个 牢房 里 的 内 犯 同时 ， 必 须要 贿赂 两 旁 相 邻 
牢房 里 的 内 犯 一 枚 金币 。 
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另外 ， 为 了 防止 释放 的 消息 在 相 邻 牢房 间 传 开 ， 不 仅 两 旁 直接 相 邻 的 
牢房 ， 所 有 可 能 听 到 消 


轧 的 内 犯 ， 即 直到 至 牢房 为 止 或 直到 监 狐 两 器 为 止 ,此 间 的 所 有 办 犯 都 
必须 给 一 枚 金币 。 


释放 后 给 金币 的 例子 


现在 要 释放 a 1, a 2, .…, aa 号 牢房 里 的 Q 名 办 犯 ， 释 放 的 顺序 还 没 确 
定 。 如 果 选 择 所 需 金 币 数 


量 尽量 少 的 顺序 释放 ， 最 少 需 要 多 少 枚 金币 ? 


限制 条 件 

Ul<N<100 

UQ<P 

Small 

[}1<P<100 

[1<0Q<5 

Large 

[] 1 <P < 10000 

[}1<Q< 100 

样 例 1 

输入 

P=8, Q=1, A= {3} 

输出 

7 〈 必 须要 给 剩 下 的 7 个 人 一 枚 金币 ) 
样 例 2 

输入 

P= 20, Q=3, A = (3, 6, 14} 
输出 


35( 按 照 牢房 14、 牢 房 6、 牢 房 3 的 顺序 释放 ， 则 需要 19 + 12 + 4E0354 
金币 ， 是 最 少 的 ) 
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DRM Ke, PNT RNAV Za, PATE oa A 
段 ， 此 后 这 两 段 束 相互 独立 了 。 1 


2 

释放 后 分 成 两 部 分 的 情况 

3 

释放 上 图 中 次 里 的 因 犯 时 

O 此 时 所 需 的 金币 数量 

O 释放 左 侧 部 分 (O) 所 需 的 金币 总 数 
4 


O 释放 右 侧 部 分 (D) 所 需 的 金币 总 数 


三 者 的 总 和 就 是 所 需 的 金币 忌 数 。 只 要 不 断 递 归 地 枚 举 最 初 释放 的 
AH LM et 


数 ， 束 能 求 出 答案 了 。 


5 


这 里 ， 递 归 计 算 过 程 中 作为 计算 对 象 的 连续 部 分 ， 其 两 端 是 空 牢房 或 
A > 因此， 作为 计算 


对 象 的 连续 部 分 一 共有 O (Q 2) 个 。 所 以 ， 利 用 动态 规划 按 顺 序 计算 ， 
就 能 够 在 O (Q 3) 时 间 内 求解 。 


6 

/ 输入 

int P, Q, AIMAX_Q + 2];/A 中 保存 输入 数据 ， 下 标 从 1 开始 
// dplillj] := 释放 (i, j) 所 需 的 金币 QQ) 


int dp[MAX_Q + 1][MAX_Q + 2]; 


7 
void solve() { 

WAT AE, RPMI AAO) 
Ar[0] = 0; 

A[Q+ 1] =P +1; 

8 

/初始 化 

for (int q = 0; q <= Q; q++) { 

dplqllq + 1] = 0; 

} 

9 


/ 从 短 的 区 间 开 始 填充 dp 

for (int w = 2; w <= Q + 1; w++) { 
for (inti=0;i+w<=Q + 1;1++) { 
/ 计算 dp[i]Dj] 

intj =it+tw,t=INT_MAX; 


10 


O dp[i][j] 表 示 的 是 ， 将 从 A 跨 号 囚犯 到 A[j] 号 因 犯 〈 不 含 两 端的 内 犯 ) 
的 连续 部 分 里 的 所 有 因 犯 都 释放 时 ,所 需 的 最 少 金 


币 总 数 。 
11 


D 为 了 更 方便 地 处 理 两 端的 情况 ， 我 们 把 左 端 当成 0 号 囚犯 ， 右 端 当成 
Q+1 写 囚犯。 这 样 ，dp[0][Q+1] 束 十 答案 。 
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/ 枚 举 最 初 释放 的 囚犯 ,计算 最 小 的 费用 
for (int k =i + 1; k <j; k++) { 


t = min(t, dplil[k] + dp[k][j]); 


} 

/ 最 初 的 释放 还 需要 与 所 释放 因 犯 无 关 的 A[j]-A[i]-2 枚 金币 
dplillj] = t + ALj]- Ali] - 2; 

} 

} 

printf("%d\n", dp[0][Q + 1]); 

} 

2.7.4 Millionaire 

Millionaire (2008 APAC local onsites C) 


你 被 邀请 到 某 个 电视 节目 中 去 玩 下 面 这 个 游戏 。 一 开始 你 有 x 元 钱 ， 
接着 进 M FORM > EE 


e, 可 以 将 所 持 的 任意 一 部 分 ) 钱 作为 赌注 。 财 注 不 光 可 以 十 整数 ， 
也 可 以 是 小 数 。 一 分 钱 


不 押 或 全 押 都 没有 关系 。 每 一 轮 都 有 的 概率 可 以 赢 ， 赢 了 赌注 束 会 
au, ¥ SEMIS T° 


如 果 你 最 后 挂 有 1000000 元 以 上 的 钱 的 话 ， 束 可 以 把 这 些 钱 种 回 家 。 
请 计算 当 你 采取 最 优 策 


略 时 ， 获 得 1000000 元 以 上 的 钱 并 带 回 家 的 概率 。 
限制 条 件 


[O<P<1.0 


D 1 < X < 1000000 


Small 


[1<M<5 

Large 

H1i<M<15 

样 例 1 

输入 

M = 1, P=0.5, X = 500000 
输出 


0.500000( 一 开始 便 全 押 ) 
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样 例 2 

1 

输入 


M = 3, P = 0.75, X = 600000 


2 

输出 
0.843750 

1. 该 问题 的 难点 


3 


续 性 ”是 这 1 个 问题 的 一 大 特点 。 每 一 轮 可 押 的 赌注 不 一 定 非 是 整 
无 限 种 可 能 ， 所 


因而 
ETI ABIL ° 


E 
a, 
以 完 


4 
2. 化 连续 为 离散 


但 是 ， 认 真 思考 一 下 束 会 发 现 ， 我 们 只 需要 检查 “无 限 种 可 能 ”中 的 “有 
限 种 可 能 * 束 足够 了 ， 


TA lear ery 下 最 后 一 轮 会 出 现 哪些 情 
Ue 


5 


D 如 果 你 持 有 1000000 元 以 上 的 钱 ， 就 没有 再 财 的 必要 了 “。 有 1 的 概率 可 
以 带 钱 回 家 。 


O 如 果 你 持 有 500000 元 以 上 的 钱 ， 不 妨 全 押 了 。 有 了 的 概率 可 以 带 钱 
回 家 。 虽 然 不 全 押 也 是 可 以 


的 ， 不 过 因为 要 求 要 达到 1000000 元 ， 所 以 不 论 怎样 都 是 这 轮 赢 了 就 能 
市 钱 回 家 ， 输 了 就 不 能 


6 
市 钱 回 家 。 


O 如 果 你 持 有 不 到 500000 元 的 钱 ， 那 已 经 无 可 奈何 了 。 有 0 的 概率 可 以 
带 钱 回 家 。 


7 


8 


9 


最 后 一 轮 带 钱 回 家 的 概率 


虽然 赌 注 是 连续 的 (有 无 限 种 可 能 ) ， 但 实际 的 概率 是 像 上 面 这 样 的 
阶梯 函数 。 因 此 ， 只 需 考虑 


10 
处 于 这 三 个 范围 中 的 哪 一 个 就 可 以 了 。 
那么 ， 最 后 两 轮 时 叉 如 何 呢 ? 同样 地 ， 这 次 也 分 为 5 种 情况 。 


11 
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最 后 两 轮 带 钱 回 家 的 概率 


某 个 范围 中 ， 即 使 所 持 的 钱 数 不 同 ， 最 后 可 以 市 钱 回 家 的 概率 也 是 完 
全 一 样 的 。 而 且 同 样 地 ，M 


轮 时 只 要 考虑 2 M +1 种 情况 束 足 够 了 。 这 样 ， 束 可 以 设计 出 穷竭 搜索 
的 算法 了 。 


3. 动态 规划 
接 下 来 ， 将 穷 疯 搜索 改 成 动态 规划 ， 就 能 够 求解 Large 了 > 
/ 输入 


int M, X; 


double P; 


double dp[2][(1 << MAX_M) + 1]; 

void solve() { 

int n = 1 << M; 

double *prv = dp[0], *nxt = dp[1]; 
memset(prv, 0, sizeof(double) * (n + 1)); 
prv[n] = 1.0; 

for (int r = 0; r < M; r++) { 

for (int i = 0; i <= n; i++) { 

int jub = min(i, n - i); 

double t = 0.0; 

for (int j = 0; j <= jub; j++) { 

t = max(t, P * prv[i + j] + (1 - P) * prv[i - jJ; 
} 

nxt[i] = t; 

} 

swap(prv, nxt); 

} 

int i = (1X * n / 1000000; 
printf("%.6f\n", prv[i]); 


} 


练习 题 135 

练习 题 

1 

2.1 REMA ARR” 

POJ 1631: Bridging signals 

U 

POJ 3666: Making the Grade 
深度 优先 搜索 

POJ 2392: Space Elevator 

POJ 1979: Red and Black 

POJ 2184: Cow Exhibition 

2 

AOJ 0118: Property Distribution 
2.4 加 工 并 存储 数据 的 数据 结构 
AOJ 0033: Ball 

POJ 3009: Curling 2.0 

0 优先 队列 

O 广度 优先 搜索 

POJ 3614: Sunscreen 


POJ 2010: Moo University - Financial Aid 


3 

AOJ 0558: Cheese 

POJ 3669: Meteor Shower 

0 并 得 集 

AOJ 0121: Seven Puzzle 

POJ 2236: Wireless Network 

O Fi wate FR 

POJ 1703: Find them, Catch them 
AOJ 2170: Marked Ancestor 

4 

POJ 2718: Smallest Difference 
POJ 3187: Backward Digit Sums 
2.5 它们 其 实 都 是 “图 ” 


POJ 3050: Hopscotch 


O 最 短路 

AOJ 0525: Osenbei 

AOJ 0189: Convenient Location 
2.2 一 往 直 前 ! 贪心 法 

5 


POJ 2139: Six Degrees of Cowvin Bacon 


D 区 间 

POJ 3259: Wormholes 

POJ 2376: Cleaning Shifts 
POJ 3268: Silver Cow Party 
POJ 1328: Radar Installation 
AOJ 2249: Road Construction 
AOJ 2200: Mr. Rito Post Office 
6 

POJ 3190: Stall Reservations 
O 其 他 

0 最 小 生成 树 

POJ 2393: Yogurt factory 
POJ 1258: Agri-Net 

POJ 1017: Packets 

POJ 2377: Bad Cowtractors 
AOJ 2224: Save your cat 

7 

POJ 3040: Allowance 

POJ 1862: Stripies 


POJ 2395: Out of Hay 


POJ 3262: Protecting the Flowers 
2.6 EM aN ] 

2.3 记 孙 结 末 再 利用 的 “动态 规划 ” 
O 轧 转 相 除 法 

8 

O 基础 的 动态 规划 算法 

AOJ 0005: GCD and LCM 

POJ 3176: Cow Bowling 

POJ 2429: GCD & LCM Inverse 
POJ 2229: Sumsets 

POJ 1930: Dead Fraction 

POJ 2385: Apple Catching 

O 素数 

9 

POJ 3616: Milking Time 

AOJ 0009: Prime Number 

POJ 3280: Cheapest Palindrome 
POJ 3126: Prime Path 

O 优化 递 推 关系 式 


POJ 3421: X-factor Chains 


POJ 1742: Coins 

POJ 3292: Semi-prime H-numbers 
10 

POJ 3046: Ant Counting 

O Rik ie Fe 

POJ 3181: Dollar Dayz 

POJ 3641: Pseudoprime numbers 
O 需 稍 加 思考 的 题目 

POJ 1995: Raising Modulo Numbers 
POJ 1065: Wooden Sticks 

11 
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3.2 常用 技巧 精 选 (一 ) 
在 此 我 们 介绍 一 些 程序 设计 竞赛 中 的 常用 技巧 。Q@ 


3.2.1 RBGEO) 
Subsequence (POJ No.3061) 


AEREA n 的 数列 整数 a 0, a 1,..., an D1 以 及 整数 S。 求 出 总 和 不 小 
于 5 的 连续 子 序列 的 长 度 的 


最 小 值 。 如 果 解 不 存在 ， 则 输出 0。 
限制 条 件 

D10<n<105 

D0<ai<104 

DS< 108 

样 例 1 

输入 


a=(5,1,3, 5, 10, 7, 4, 9, 2, 8) 
输出 

2 (5+10) 

样 例 2 

输入 


n=5 


@ 本 世 所 涉及 的 多 种 技巧 ， 其 应 用 范围 都 比 这 里 所 介绍 的 模型 要 广 得 
多 ， 注 意 强 含 在 技巧 内 的 算法 思想 。 许 多 看 似 非 沼 


复 杀 困难 的 问题 ， 其 问题 天 键 都 可 以 运用 这 里 看 似 商 单 的 技巧 处 理 。 


译 者 注 


@ 这 里 直接 使 用 了 日 文 原文 的 汉字 ， 斥 取 法 通 音 羡 指 对 数组 保存 一 对 
下 标 (起 点 、 终 点) ， 然 后 根据 实际 情况 交替 推进 


两 个 端点 直到 得 出 答案 的 方法 ， 这 种 操作 很 像 是 尺 凤 (日 文中 称 为 尺 
WE) 扑 行 的 方式 故 得 名 。 一 一 译 者 注 
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S=11 


1 
a={1,2, 3,4,5} 
输出 

3 (3+4+5) 

2 


由 于 所 有 的 元 素 都 大 于 零 ， 如 果子 序列 [ s ,t) 满 足 as+...+atDl>S ， 
那么 对 于 任何 的 t<t 一 定 有 as +... 


+ at'D1>S。 此 外 对 于 区 间 [s,t) 上 的 总 和 来 说 如 末 令 


3 


sum(i )=a0+a1+...+ ai[]1 
ALA 


as + as +1+...+ at []1=sum( t )[|sum(s ) 4 


因此 预先 以 O(n ) 的 时 间 计 算 好 sum 的 话 ， 就 可 以 以 O (1) 的 时 间 计 算 
区 间 上 的 总 和 。 这 样 一 来 ， 子 


序列 的 起 点 s 确定 以 后 ， 便 可 以 用 二 分 搜索 快速 地 确定 使 序列 和 不 小 
F S 的 结尾 上 的 最 小 值 。 


5 

/输入 

int n, S; 

int alMAX_N]; 

int sum[MAX_N + 1]; 

6 

void solve() { 

/ 计算 sum 

for (int i = 0; i < n; i++) { 
sum[i + 1] = sum[i] + ali]; 
7 

} 

if (sum[n] < S) { 

/ 解 不 存在 
printf("0\n"); 

8 


return; 


} 

int res = n; 

for (int s = 0; sum[s] + S <= sum[n]; s++) { 
9 

/利用 二 分 搜索 求 出 t 


int t = lower_bound(sum + s, sum + n, sum[s] + S) - sum; res = min(res, t - 


S); 

} 

printf("%d\n", res); 
10 

} 


这 个 算法 的 复杂 度 是 O (nlogm)， 虽 然 足以 解决 这 个 问题 ， 但 我 们 还 
可 以 更 加 高 效 地 求解 。 我 们 设 


以 as 开始 总 和 最 初 大 于 S 时 的 连续 子 序列 为 as +...+ at 
11 


01, 325 


315 0 7419| 2 
\ 
spapafsjopopejofojo 
\ o 
了 5135 和 7 492 
\ 
BRED 
NN 
spapafsfoprfojolofo 
Se 
spapapspoprjefojojo 
N N 


L 1 2 MC | ant 7 1 0 7 Q 
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as +1+...+ at [|2<as+...+at[l2<5S 


所 以 从 as +1 开 始 总 和 最 初 超过 S 的 连续 子 序列 如 果 是 as +1+...+ at TI 
的 话 ， 则 必然 有 tst。 利 用 这 一 性 


质 便 可 以 设计 出 如 下 算法 : 
(1) U s = t = sum = 0 初始化。 
(2) 只 要 依然 有 sum< S ， 束 不 断 将 sum 增加 at ， 并 将 上 增加 1。 


(3) 如 采 (2) 中 无 法 满足 sum> S 则 终止 。 否 则 的 话 ， 更 新 res = min(res, t 
Us)。 


(4) 将 sum RE as, ，s 增加 1 然后 回 到 (2) 。 


对 于 这 个 算法 ， 因 为 1 最 多 变化 n 次 ， 因 此 只 需 O(n ) 的 复 洒 度 束 可 以 
求解 这 个 问题 了 。 


void solve() { 


intres=n+1; 

int s = 0, t= 0, sum = 0; 
for (53) { 

while (t < n && sum < S) { 
sum += alt++J; 

} 

if (sum < S) break; 

res = min(res, t - s); 


sum -= a[s++]; 


} 

if (res > n) { 

/ 解 不 存在 

res = 0; 

} 

printf("%d\n", res); 
} 


样 例 数据 1 对 应 的 区 间 的 变化 


3.2 常用 技巧 精 选 (一 ) 149 


像 这 样 反 复 地 推进 区 间 的 开头 和 末尾 ， 来 求 取 满 足 条 件 的 最 小 区 间 的 
方法 被 称 为 尺 取 法 。 


1 


Jessica’s Reading Problem (POJ No.3320) 


为 了 准备 考试 ，Jessica 开始 读 一 本 很 厚 的 读本 。 要 想 通 过 考试 ， 必 须 
把 课本 中 所 有 的 知识 所 


2 


都 掌握 。 这 本 书 总 共有 PP 页， 第 i 页 恰好 有 一 个 知识 点 qi (每 个 知识 
点 都 有 一 个 整数 编号 ) 。 


ee a 点 可 能 会 被 多 次 提 到 ， 所 以 她 希望 通过 阅读 其 中 连 
续 的 一 些 页 把 所 有 的 知 


识 点 都 覆盖 到 。 给 定 每 页 写 到 的 知识 点 ， 请 求 出 要 阅读 的 最 少 页 数 。 
3 
限制 条 件 


[]1<P< 106 


4 
样 例 


a 

输出 

2 (只 要 阅读 第 1 页 和 第 2 页 即 可 ) 
6 


我 们 假设 从 某 一 页 s 开始 阅读 ， 为 了 覆盖 所 有 的 知识 点 需要 阅读 到 t。 
这 样 的 话 可 以 知道 如 采 从 s +1 


那么 必须 阅读 到 t'>t 页 为 止 。 由 此 这 题 也 可 以 使 用 尺 
IE © 


7 


在 某 个 区 间 [s ,t] 已 经 窗 盖 了 所 有 的 知识 后 的 情况 下 ， 下 一 个 区 间 [s 
+1, 4"](t'>t) 要 如 何 求 出 呢 ? 


所 有 的 知识 点 都 被 履 兰 [每 个 知识 点 出 现 的 次 数 都 不 小 于 1 
8 


由 以 上 的 等 价 和 关系， 我 们 可 以 用 二 又 树 等 数据 结构 来 存储 [s e] IAE 
每 个 知识 点 的 出 现 次 效 ， 这 


样 把 最 开头 的 页 s 去 除 后 便 可 以 判断 [ s +1, t ETER ° 


从 区 间 的 最 开头 把 s 取 出 之 后 ， 页 s 上 书写 的 知识 点 的 出 现 次 数 就 要 减 
一 ， 如 采 此 时 这 个 知识 点 的 


9 


出 现 次数 为 0 了 ， 在 同一 个 知识 点 再 次 出 现 前 ， 不 停 将 区 间 末 尾 上 问 后 
推进 即 可 。 每 次 在 区 间 林 尾 


人 退 加 页 t 时 将 页 t 上 的 知 误 后 的 出 现 次 数 加 1， 这 样 束 完成 了 下 一 个 区 
间 上 各 个 知识 点 出 现 次数 的 更 


新 ， 通 过 重复 这 一 操作 可 以 以 O (PlogP ) 的 复杂 度 求 出 最 小 的 区 间 。 
10 


/输入 


int P; 

int a MAX P]; 

void solve() { 

11 

/ 计算 全 部 知识 点 的 总 数 
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set<int> all; 
for (int i = 0; i < P; i++) { 
all.insert(a[i]); 


} 


int n = all.size(); 

/利用 扩 取 法 来 求解 

int s = 0, t= 0, num = 0; 
map<int, int> count; // 知识 点 = 出 现 次 数 的 映射 
int res = P; 

for (;;) { 

while (t < P && num <n) { 

if (count[a[t++]]++ == 0) { 

/ 出 现 新 的 知识 点 

num++; 

} 

} 

if (num < n) break; 

res = min(res, t - s); 

if (--count[a[s++]] == 0) { 

1/ 某 个 知识 护 的 出 现 次 数 为 0 了 
num--; 

} 

} 


printf("%d\n", res); 


} 
3.2.2 反 转 (开关 问题 ) 
Face The Right Way (POJ No. 3276) 


N AEIR T + BA FA AA a o A SULA AaB E 
BIT, KRAK T 


—8 BIRMANIA LIA MUAK, N 
BRE UK te (EK 


头 连 续 的 牛 转向 。 请 求 出 为 了 让 所 有 的 牛 都 能 面向 前 方 需要 的 最 少 的 
操作 次 数 M 和 对 应 的 


最 小 的 K。 

限制 条 件 

[] 1< N <5000 

样 例 

输入 

N=7 

BBFBFBB F: 面向 前 方 ，B: 面向 后 方 ) 


aele 
Heeel 
Hehee 


wo 
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输出 


1 


M=3 
( 先 反 转 1~3 号 的 三 头 牛 ， 然 后 再 反 转 3~5 号 ， 最 后 反 转 5~7 号 ) 


3 
4 
K=3 时 的 解法 


5 


首先 我 们 来 看 看 对 于 一 个 特定 的 K 如 何 求 出 让 所 有 的 牛 面 朝 前 方 的 最 
小 操作 次 数 。 如 来 把 牛 的 方 


回 作 为 状态 进行 搜索 的 话 ， 由 于 状态 数 有 2 N 个 ， 古 无 法 在 时 限 内 找 出 
答案 的 。 那 么 不 搜索 的 话 要 


6 
EA DIE? 


首先 ， 交换 区 间 反 转 的 顺序 对 结果 是 没有 影响 的 。 此 外 ， 可 以 知道 对 
同一 个 区 间 进 行 两 次 以 上 


的 反 转 是 多 余 的 。 由 此 ， 问 题 束 转化 成 了 求 需要 被 反 转 的 区 间 的 集 
5° TERMAS DN 


7 


左 端的 牛 。 包 含 这 头 牛 的 区 间 只 有 一 个 ， 因 此 如 果 这 头 牛 面 朝 前 方 ， 
我 们 就 能 知道 这 个 区 间 不 


需要 反 转 。 


反之 ， 如 采 这 头 牛 面 瑚 后方， 对 应 的 区 间 吏 必须 进行 反 转 了 “。 而 且 在 
此 之 后 这 个 最 左 的 区 间 就 再 


8 


也 不 需要 考虑 了 。 这 样 一 来 ， 通 过 首先 考虑 最 左 端的 牛 ， 问 题 的 规模 
忠 缩 小 了 人 1。 不 断 地 重复 下 


去 ， 就 可 以 无 需 搜 索 求 出 最 少 所 需 的 反 转 次 数 了 。 


此 外 ， 通 过 上 面 的 分 析 可 以 知道 ， 名 略 掉 对 同一 个 区 间 重 复 反 转 这 类 
多 余 操 a 作 之 后 ， 只 要 存在 让 


9 
所 有 的 牛 都 朝 前 的 方法 ， 那 么 操作 就 和 顺序 无 关 可 以 唯一 确定 了 。 


这 个 算法 的 复杂 度 又 如 何 呢 ? 首先 我 们 需要 对 所 有 的 天 都 求解 一 次 ， 
对 于 每 个 天 我 们 都 要 从 最 左 


端 开 始 来 考虑 N 头 牛 的 情况 。 此 时 最 坏 情况 下 需要 进行 WDK+1 次 的 
反 转 操作 ， 而 每 次 操作 又 要 反 


10 


转 天 头 牛 ， 于 是 总 的 复杂 度 就 是 O ( N3)。 这 样 的 话 还 不 足以 在 时 限 内 
解决 问题 。 但 是 区 间 反 转 的 部 


分 还 是 很 容易 进行 优化 的 。 

fi := 区间 [i,i+K0D1] 进 行 了 反 转 的 话 则 为 1， 否 则 为 0 
11 

152 第 3 Æ 出 类 拔 茎 一 一 中 级 篇 

这 样 ， 在 考虑 第 头 牛 时 ， 如 果 


¡1 

U 

0 

Ali] 为 奇数 的 话 ， 则 这 头 牛 的 方 癌 与 起 始 方 网 是 相反 的 ， 
JUiUK1 

0 


否则 方向 不 变 。 由 于 
i 
¡1 


U 


UflsiIOUF A OFLiOFLIOK 01) 
iDi 


iDipK1 
U 


所 以 这 个 和 每 一 次 都 可 以 用 常数 时 间 计 算出 来 ， 复杂 度 就 降 为 了 O(N 
2)， 能 在 时 限 内 解决 了 。 


1/ 输入 

int N; 

int dir[MAX_N]; // 牛 的 方向 (0:E 1:B) 

int f[MAX_N]; // 区 间 [Li+K-1] 是 人 否 进行 反 转 
/固定 K， 求 对 应 的 最 小 操作 回 数 

/无 解 的 话 则 返回 -1 

int calc(int K) { 

memset(f, O, sizeof(f)); 

int res = 0; 


int sum = 0; / fA #0 


for (int i = 0; i + K <= N; i++) { 

/ 计算 区 间 [i,it+K-1] 

if ((dir[i] + sum) % 2 != 0) { 
AT 

res++; 

f[i] = 1; 

} 

sum += f[i]; 

if(i-K+1>=0){ 

sum -= f[i- K + 1]; 

} 

} 

"RER NUERA TE SB EN BY DL, 
for (inti=N-K+1;i<N;i++) { 
if ((dir[i] + sum) % 2 != 0) { 
/无 解 

return -1; 

} 

if (i- K+1>=0) { 


sum -= f[i- K + 1]; 


} 

} 

return res; 

} 

void solve() { 
intK=1,M=N; 


for (int k = 1; k <= N; k++) { 


int m = calc(k); 
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if (m >= 0 && M > m) { 


1 

=m; 
K=k; 
} 
} 


printf("%d %d\n", K, M); 


} 
2 
Fliptile (POJ No.3279) ® 
3 


农夫 约翰 知道 聪明 的 牛 产 奶 多 。 于 是 为 了 提高 牛 的 智商 他 准备 了 如 下 
游戏 。 有 一 个 MxNN 的 


格子 ， 每 个 格子 可 以 翻转 正 反 面 ， 它 们 一 面 古 黑色 ， 男 一 面 是 日 色 。 
REITE TE a DE 


日 色 ， 日 色 的 格子 翻 罗 过 来 则 是 墨色。 游戏 要 做 的 吏 是 把 所 有 的 格子 
都 翻转 成 白色 。 不 过 因 


4 


为 牛 蹄 很 大 ， 所 以 每 次 翻转 一 个 格子 时 ， 与 它 上 下 左右 相 邻 接 的 格子 
BRA > AA RIS 


子 太 麻烦 了 ， 所 以 牛 都 想 通 过 尽 可 能 少 的 次 数 把 所 有 格子 都 翻 成 日 
色 。 现 在 给 定 了 每 个 格子 


的 颜色 ， 请 求 出 用 最 小 步 数 完成 时 每 个 格子 翻转 的 次 数 。 最 小 步 数 的 
解 有 多 个 时 ， 输 出 字典 


序 最 小 的 一 组 。 解 不 存在 的 话 ， 则 输出 IMPOSSIBLE > 
5 
限制 条 件 


|1<M,N<15 


6 
样 例 


输入 

7 

M=4 

N=4 

每 个 格子 的 颜色 如 下 (0 表示 白色 ，1 表 示 黑 色 ) 


1001 


OIT BAF A AT Tal, BERN TExtended 
Lights Out (Greater New York 2002) 这 道 题 。 高 斯 消 


元 法 也 可 以 用 于 求解 一 组 可 行 解 ， 并 且 通 过 这 些 分 析 可 以 知道 ， 目 由 
变 元 的 个 数 不 会 超过 NN 个 ， 所 以 也 可 以 用 于 求解 


11 
最 优 解 。 一 一 译 者 注 


1000 IN 


EE Fy De 


解法 示意 


目 完 ， 同 一 个 格子 翻转 两 次 的 话 束 会 恢复 原状 ， 所 以 多 次 翻转 是 多 余 
的 。 此 外 ， 翻 转 的 格子 的 集 


合 相同 的 话 ， 其 次 序 是 无 关 紧 要 的 。 因 此 ， 总 共有 2 NM 种 翻转 的 方 
法 。 不 过 这 个 解 空 间 太 大 了 ， 


我 们 需要 想 出 更 有 效 的 办 法 。 


让 我 们 再 回顾 一 下 前 面 的 问题 。 在 那 道 厦 中 ， 让 最 左 闹 的 牛 反 转 的 方 
法 只 有 1 种 ， 于 是 用 直接 判 


断 的 方法 确定 惑 可 以 了 。 同 样 的 方法 在 这 里 还 行 得 通 吗 ? 


不 妨 先 看 看 最 左上 角 的 格子 。 在 这 里 ， 除 了 翻转 (11) 之 外 ， 翻 转 (1,2) 和 
(2,1) 也 可 以 把 这 个 格子 翻 


转 ， 所 以 像 之 前 那样 直接 确定 的 办 法 行 不 通 。 


于 走 不 妨 先 指定 好 最 上 面 一 行 的 翻 胃 方法 。 此 时 能 够 翻 力 (1.1) 的 只 剩 
下 (2,1) 了 ， 所 以 可 以 直接 判 


MODERA ZE o MIO DAA, N ) 都 能 这 样 判断 ， 如 此 反复 下 
去 就 可 以 确定 所 有 格子 的 翻 


转 方 法 。 最 后 ( M ,1)~( M,N ) 如 采 并 非 全 为 日 色 ， 束 意味 着 不 存在 可 行 
的 操作 方法 。 


像 这 样 ， 先 确定 第 一 行 的 翻转 方式 ， 然 后 可 以 很 容易 判断 这 样 是 否 存 
在 解 以 及 解 的 最 小 步 数 十 多 


少 ， 这 桂 将 第 一 行 的 所 有 翻转 方式 都 笠 试 一 次 束 能 求 出 整个 问题 的 最 
小 步 数 。 这 个 算法 中 最 上 面 


一 行 的 翻转 方式 共有 2 N 种 ， 复 杂 度 为 O (MN2N)。 
/ 邻接 的 格子 的 坐标 


const int dx[5] = {-1, 0, 0, 0, 1}; 


const int dy[5] = {0, -1, 0, 1, O}; 

// 输入 

int M, N; 

int tile[MAX_MI[MAX_N]; 

int opt[MAX_M][MAX_N]; // 保存 最 优 解 
int flip[MAX_M][MAX_N]; // 保存 中 间 结 果 
/ BA (xy WE 

int get(int x, int y) { 

int c = tile[x][y]; 

for (int d = 0; d < 5; d++) { 

int x2 = x + dx[d], y2 = y + dy[d]; 

if (0 <= x2 && x2 < M && 0 <= y2 && y2 < N) { 
c += flip[x2][y2]; 

} 

} 

return c % 2; 
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} 

1 

/ 求 出 第 1 行 确定 情况 下 的 最 小 操作 次 数 


/ 不 存在 解 的 话 返回 -1 

int calc() { 

/ 求 出 从 第 2 行 开 始 的 翻转 方法 
for (int i = 1; i < M; i++) { 

2 

for (int j = 0; j < N; j++) { 


if (get(i - 1, j) != 0) { 


/ -1Lj) 是 黑色 的 话 ， 则 必须 翻转 这 


位 ip[iU] = 1; 

} 

3 

h 

} 

Ft Bate eS A 
for (int j = 0; j < N; j++) { 
4 

if (get(M - 1, j) !=0) 4 
/无 解 

return -1; 


} 


个 格子 


} 

5 

// 统计 翻转 的 次 数 

int res = 0; 

for (int i = 0; i < M; i++) { 
for (int j = 0; j < N; j++) { 
6 

res += fliplillj]; 

} 

} 

return res; 

} 

7 

void solve() { 

int res = -1; 
PREM I — ITA PE 


for (int i = 0; i < 1 << N; i++) { 


8 
memset(flip, 0, sizeof(flip)); 


for (int j = 0; j < N; j++) { 


flip[olIN -j - 1] =i>>j &1; 

} 

int num = calc(); 

9 

if (num >= 0 && (res < 0 || res > num)) { 
res = num; 

memcpy(opt, flip, sizeof(flip)); 

} 

} 

10 

if (res < 0) { 

/无 解 

printf("IMPOSSIBLE\n"); 

} else { 

11 

for (int i = 0; i < M; i++) { 
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for (intj = 0; j < N; j++) { 


printf("%d%c", optlilljl, j + 1 == N ?'\n': ''); 


} 


} 
} 
} 
专 位 集合 的 整数 表示 


前 面 的 代码 里 ， 为 了 笑 试 第 一 行 的 所 有 可 能 性 ， 使 用 了 集合 的 整数 表 
现 。 在 程序 中 表示 集合 


的 方法 有 很 多 种 ， 当 元 素数 比较 少时 ， 像 这 样 用 二 进 制 码 来 表示 比较 
方便 。 集 合 {0,1...,nD1} 


的 子 集 S 可 以 用 如 下 方式 编码 成 整数 。 


f(S)O02i 

is 

O 

像 这 样 表示 之 后 ， 一 些 集合 运算 可 以 对 应 地 写成 如 下 方式 。 
BN ner er mtr eax 0 
DRERR TTL a EN ann 1<<i 
DEEE n 个 元 素 的 集合 {0,1,..., n Gb: ec. (1<<n)-1 

0 判断 第 i 个 元 素 是 否 属于 集合 S :ee if (S>>i&1) 
0 向 集合 中 加 入 第 i 个 元 素 S { 

OD euere SI1<<i 

E ZH NEE SIT ee S&~(1<<i) 


DRE SFA TAFSER SU Ts iseina SIT 


a ai S&T 


此 外 ， 想 要 将 集合 {0,1,.…., n 01} 的 所 有 子 集 枚 举 出 来 的 话 ， 可 以 像 下 面 
这 样 书写 


for (int S = 0; S < 1 << n; S++) { 
// 对 于 集 的 处 理 
} 


按照 这 个 顺序 进行 循环 的 话 ， 5 便 会 从 空 集 开始 ， 然 后 按照 {0}、 
10,1. Oat nh 0} 


的 升序 顺序 枚 举 出 来 。 


接 下 来 介绍 一 下 如 何 枚 举 某 个 集合 sup 的 子 集 。 这 里 sup 是 一 个 二 进 制 
码 ， 其 本 映 也 是 某 


个 集合 的 子 集 。 例 如 给 定 了 01101101 这 样 的 集合 ， 要 将 01100000 或 
者 00101101 等 子 集 枚 举 


出 来 。 前 面 是 从 0 开始 不 断 加 1 来 枚 举 出 了 全 部 的 子 集 。 此 时 ，sub+1 
并 不 一 定 是 sup 的 


子 集 。 而 (sub+1)&sup 虽然 是 sup 的 子 集 ， 可 是 很 有 可 能 依旧 是 sub, 
没有 任何 改变 。 


所 以 我 们 要 反 过 来 ， 从 sup 开始 每 次 减 1 直到 0 为止 由 于 sub-1 并 不 


一 定 是 sup 的 子 集 ， 


所 以 我 们 把 它 与 sup 进行 按 位 与 &。 这 样 的 话 束 可 以 将 sup 所 有 的 子 集 
按照 降序 列举 出 来 。 


(sub-1)&sup ZA sup 中 的 0 而 从 sub 中 减 去 1。 
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int sub = sup; 


1 
do { 
// 对 子 集 的 处 理 


sub = (sub - 1) & sup; 


} while(sub != sup); // 处 理 完 0 之 后 ， 会 有 -1&sup=sup 
2 


最 后 我 们 介绍 一 下 枚 举 {0,1,.…, n D1} 所 包含 的 所 有 大 小 为 k 的 子 集 的 方 
法 。 通 过 使 用 位 运算 我 


们 可 以 像 如 下 代码 所 示人 简单 地 按照 字典 序 升 序 地 枚 举 出 所 有 满足 条 件 
的 二 进 制 码 。 


int comb = (1 << k) - 1; 


while (comb < 1 <<n) { 

3 

/这 里 进行 针对 组 合 的 处 理 

int x = comb & -comb, y = comb + x; 
comb = ((comb & ~y)/ x >> 1) | y; 

} 

4 


按照 字典 序 的 话 ， 最 小 的 子 集 是 (1<<k)-1， 所 以 用 它 作 为 初始 值 。 现 在 
我 们 求 出 comb 其 


后 的 二 进 制 码 。 例 如 0101110 之 后 的 是 0110011，0111110 之 后 的 是 
1001111。 下 面 是 求 出 comb 下 一 个 二 进 制 码 的 方法 。 


5 
(1) 求 出 最 低位 的 1 开始 的 连续 的 1 的 区 间 (0101110 > 0001110) 


(2) 将 这 一 区 间 全 部 变 为 0， 并 将 区 间 左 侧 的 那个 0 变 为 1 
(0101110 > 0110000) 


(3) 将 第 (1) 步 里 取出 的 区 间 右 移 ， 直 到 剩 下 的 1 的 个 数 减 少 了 1 个 
(0001110 > 0000011) 


(4) 将 第 (2) 步 和 第 (3) 步 的 结果 按 位 取 或 (0110000|0000011=0110011) 
6 


对 于 非 零 的 整数 ，x&(-x) 的 值 束 是 将 其 最 低位 的 1 独立 出 来 后 的 值 。 这 
是 由 于 计算 机 中 人 负 


a A 
1 o 


7 


x 
x 的 二 进 制 
Ox 的 二 进 制 
x &- x 

1 0001 

1111 

0001 

2 0010 


1110 


6 0110 
1010 
0010 
70111 
1001 
0001 


将 最 低位 的 1 取出 后 ， 设 它 为 x。 那 么 通过 计算 y=comb+x, BH comb 
从 最 低位 的 1 开始 


9 


的 连续 的 1 都 置 0 了。 我 们 来 比较 一 下 ~y 和 comb。 在 comb 中 加 上 x 
后 没有 变化 的 位 ， 在 


y 中 全 都 取 相 反 的 值 。 而 最 低位 1 开始 的 连续 区 间 在 ~y 中 依然 是 1， 
区 间 左 侧 的 那个 0 在 


~y 中 也 依然 是 0。 于 是 通过 计算 z=comb&~y 就 得 到 了 最 低位 1 开始 的 
连续 区 间 。 比 如 如 采 


10 


comb=0101110， 则 x =0000010，y =0110000，z =0001110。 


同时 ，y 也 恰好 二 第 (2) 步 要 求 的 值 。 那 么 首先 将 z 不 断 石 移 ， 直 到 最 低 
位 为 1， 这 通过 计算 


z/x 即 可 完成 。 这 样 再 将 z/x 右 移 1 位 就 得 到 了 第 (3) 步 要 求 的 值 。 这 样 
BANDOS T comb 


11 
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之 后 的 下 一 个 二 进 制 列 。 因 为 是 从 n 个 元 素 的 集合 中 进行 选择 ， 所 以 
comb 的 值 不 能 大 于 等 


于 1<<n。 如 此 一 来 ， 就 完成 了 大 小 为 k 的 所 有 子 集 的 枚 举 。 


除 上 述 例子 之 外 ， 还 可 以 利用 位 运算 完成 满足 其 他 条 件 的 集合 的 枚 
举 ， 例 如 不 包含 相 邻 元 素 


的 集合 等 。 

3.2.3 弹性 碰撞 

Physics Experiment (POJ No.3684) 

用 NN 个 半径 为 R 厘 米 的 球 进行 如 下 实验 。 


在 五 米 高 的 位 置 设置 了 一 个 圆 简 ， 将 球 垂 直 放 入 (从 下 向 上 数 第 i 个 
球 的 确 端 距离 地 面 高 度 


NH+2R) 。 实 验 开 始 时 最 下 面 的 球 开 始 掉 落 ， 此 后 每 一 秒 又 有 一 个 
A Mir ETRE 


I, HRE EREHE 1 AY alt ee E EREE 。 


请 求 出 实验 开始 后 工 秒 种 时 每 个 球 底 端 的 高 度 。 假 设 重力 加 速度 为 g 
=10m/s2 ° 


限制 条 件 


[}1<N< 100 


[ 1< H < 10000 
[}1<R< 100 

[] 1 <T<10000 
样 例 1 

输入 

N=1 

H = 10 

R=10 


T = 100 


输出 
4.95 
样 例 2 
输入 


N=2 


T = 100 
输出 

2 

4.95 10.20 
3 

4 

5 

初始 状态 


首先 考虑 一 下 只 有 一 个 球 的 情形 。 这 时 只 是 单纯 的 物理 问题 。 从 高 为 
电 的 位 置 下 落 的 话 需 要 化 费 


6 
的 时 间 征 
2H 

au 

7 


g 
这 样 的 话 ， 在 时 刻 T 时 ， 令 为 满足 kt <T 的 最 大 整数 ， 那 么 
8 
9 


授 下 来 再 考虑 多 个 球 的 情形 。 乍 一 看 ， 因 为 多 个 球 之 间 会 有 碰撞 ， 必 
须 对 物理 运动 进行 模拟 ， 事 


实 上 并 没有 这 个 必要 。 我 们 来 回忆 一 下 此 前 热 号 题目 “Ants”。 在 那 道 题 
目 中 两 只 蚂蚁 相遇 后 并 


DER MERA MARE RA PEA A A 


10 


XENA CAA AIA > CICS — PR =0 的 情况 。 如 
果 认 为 所 有 的 球 都 是 一 样 的 ， 


束 可 以 无 视 它们 的 碰撞 ， 视 为 直接 互相 穿 过 继续 运动 。 由 于 在 有 磁 撞 
时 球 的 顺序 不 会 发 生 改变 ， 


11 


所 以 忽略 碰撞 ， 将 计算 得 到 的 坐标 进行 排序 后 ， 就 能 知道 每 个 球 的 最 
终 位 置 。 那 么 ，R >0 时 要 起 
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ZL Ope? 这 种 情况 下 的 处 理 方 法 基本 相同 ， 对 于 从 下 方 开 始 的 第 i 个 
球 ， 在 按照 R =0 计 算 的 结果 上 


加 上 2 Ri 就 好 了 。 
const double g = 10.0; / 重力 加 速度 
1/ 输入 


int N, H, R, T; 


double y[MAX_N]; / 球 的 最 终 位 置 
/ 求 出 T 时 刻 球 的 位 置 


double calc(int T) { 


if (T < 0) return H; 
double t = sqrt(2 * H / g); 
int k = (int)(T / t); 

if (k % 2 == 0) { 
double d= T -k * t; 
return H-g*d*d/2; 

} else { 
double d =k *t +t-T; 
return H -g*d*d/2; 

} 

} 

void solve() { 

for (int i = 0; i < N; i++) { 
yLi] = calc(T - i); 

} 

sort(y, y + N); 


for (int i = 0; i < N; i++) { 


printt("%.2f%c", yli] + 2 * R *i/100.0,i+ 1==N?'\n':"'); 
} 

} 

3.2.4 折 半 枚 举 (双向 搜索 ) O 

4 Values whose Sum is 0 (POJ No.2785) 


给 定 各 有 个 整数 的 四 个 数列 A、B、C、D。 要 从 每 个 数列 中 各 取 
出 工 个 数 ， 使 四 个 数 的 和 为 


0。 求 出 这 样 的 组 合 的 个 数 。 当 一 个 数列 中 有 多 个 相同 的 数字 时 ， 把 它 
们 作为 不 同 的 数字 看 待 。 


限制 条 件 
[1<n<4000 
O |( 数 字 的 值 | < 228 


O 本 廊 所 介绍 的 折 半 枚 举 与 传统 的 双向 搜索 并 不 相同 ， 但 其 思想 来 源 
于 传统 的 双 同 搜索 ， 有 时 候 我 们 也 会 用 双 同 搜索 来 


指 代 它 ， 特 此 注 明 。 一 一 译 者 注 
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样 例 

1 

输入 

n=6 


A = {-45, -41, -36, -36, 26, -32} 


B = {22, -27, 53, 30, -38, -54} 

2 

C = {42, 56, -37, 275,410, -6} 

D = {-16, 30, 77, -46, 62, 45} 

输出 

3 

5(-45-27+42+30=0, 26+30-10-46=0, -32+22+56-46=0, -32+30-75+77=0, 
-32-54+56+30=0) 这 个 问题 是 热 喘 题目 “抽签 ”的 复习 。 从 4 个 数列 中 选择 
的 话 总 共有 n 4 种 情况 ， 所 以 全 都 判断 一 


遍 不 可 行 。 不 过 将 它们 对 半分 成 AB 和 CD 再 考虑 ， 就 可 以 快速 解决 
了 。 从 2 个 数列 中 选择 的 话 只 有 


4 


n 2 种 组 合 ， 所 以 可 以 进行 枚 举 。 先 从 A、B 中 取出 a、b 后 ,为 了 使 
总 和 为 0 则 需要 从 C、DD 中 取出 


c+d-Jalb ° AEKA C > D 中 取 数 字 的 n2 种 方法 全 都 枚 举 出 
来 ， 将 这 些 和 排 好 序 ， 这 样 就 可 以 


运用 二 分 搜索 了 。 这 个 算法 的 复杂 度 是 O (n 2logn)。 
5 


有 时 候 ， 问 题 的 规模 较 大 ， 无 法 枚 举 所 有 元 素 的 组 合 ， 但 能 够 枚 举 一 
半 元 素 的 组 合 。 此 时 ， 将 问 


题 拆 成 两 半 后 分 别 枚 举 ， 再 合并 它们 的 结果 这 一 方法 往往 非常 有 效 。 
/ 输入 


6 


int n; 


int AIMAX_N], BIMAX_N], CIMAX_N], DIMAX_N]; 


int CD[MAX_N * MAX_N]; /C 和 D 中 数字 的 组 合 方法 
void solve() { 

7 

/ 枚 举 从 C 和 D 中 取出 数字 的 所 有 方法 
for (int i = 0; i < n; i++) { 

for (int j = 0; j < n; j++) { 

CD[i * n + j] = Cli] + Dj]; 

} 

8 

h 

sort(CD, CD+n*n); 

long long res = 0; 

for (int i = 0; i < n; i++) { 

9 

for (int j = 0; j < n; j++) { 

int cd = -(Alil + BIj]); 

/取出 C 和 D 中 和 为 cd 的 部 分 


res += upper_bound(CD, CD + n * n, cd) - lower_bound(CD, CD +n * n, 
cd); 


} 

10 

} 

printf("%lld\n", res); 
} 

11 
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超大 背包 问题 


有 重量 和 价值 分 别 为 wi, vi 的 n 个 物品 。 从 这 些 物 品 中 挑选 总 重量 不 
超过 W 的 物品 ， 求 所 有 


挑选 方案 中 价值 总 和 的 最 大 值 。 
限制 条 件 


[1<n<40 


[] 1 <wi, vi < 1015 
[]1<W< 1015 
样 例 

输入 


n=4 
w=12,1:8,2) 
v=18,2,4:2] 

W=5 

输出 

7( 挑 选 0、1、3 号 物品 ) 


这 个 问题 是 第 二 半 介 绍 过 的 背包 问题 。 不 过 这 次 价值 和 重量 都 可 以 是 
非常 大 的 数值 ， 相 比 之 下 n 


比较 小 。 使 用 DP 求解 背包 问题 的 复 洒 度 是 O(nW )， 因 此 不 能 用 来 解 
决 这 里 的 问题 。 此 时 我 们 应 该 


利用 n 比较 小 的 特点 来 寻找 其 他 办 法 。 


挑选 物品 的 方法 总 共有 2 n 种 ， 所 以 不 能 直接 枚 举 ， 但 是 像 前 面 一 样 拆 
成 两 半 之 后 再 枚 举 的 话 ， 


为 每 部 分 只 有 20 个 所 以 是 可 行 的 。 利 用 拆 成 两 半 后 的 两 部 分 的 价值 和 
重量 ， 我 们 能 求 出 原先 的 问 


题 吗 ? 我 们 把 前 半 部 分 中 的 选取 方法 对 应 的 重量 和 价值 总 和 记 为 w1、 
v1。 这 样 在 后 半 部 分 寻找 总 


重 w2< WDw1 时 使 v2 最 大 的 选取 方法 就 好 了 。 


因此 ， 我 们 要 思考 从 枚 举 得 到 的 ( w 2, v 2) 的 集合 中 高 效 寻 找 max{ v 2] w 
2< WHITE > Hd, E 

然 我 们 可 以 排除 所 有 w2[i]<w2[j] 并 且 v2[i]>v2[j] 的 j。 这 一 点 可 
以 按照 w2、v 2 的 字典 序 排序 后 


简单 做 到 。 此 后 简 余 的 元 素 都 满足 w2[i]<w2[j]v2[i]<v2[j]， 要 
计算 max{ v 2] w2< W'} 的 话 ， 只 


要 寻找 满足 w2[i]sx 的 最 大 的 i 就 可 以 了 。 这 可 以 用 二 分 搜索 完 
成 ， 剩 余 的 元 素 个 数 为 M 的 话 ， 


一 次 搜索 需要 O (log M ) 的 时 间 。 因 为 M <2(n/2)， 所 以 这 个 算法 总 的 
复杂 度 是 O (2(n /2)n)， 可 以 在 时 


限 内 解决 这 个 问题 。 

3.2 常用 技巧 精 选 (一 ) 163 
typedef long long II; 

1 

/输入 

int n; 

Il wWIMAX_N], v[MAX_N]; 

ll W; 

2 

pair<ll, 11> ps[1 << (MAX_N/ 2)]; // (重量 , 价值 ) void solve() { 
/ 枚 举 前 半 部 分 
intn2=n/2; 

3 

for (int i = 0; i < 1 << n2; i++) { 
ll sw = 0, sv = 0; 

for (int j = 0; j < n2; j++) { 


if(i>j81)( 


sw += W[j]; 

4 

sv += vlj); 

} 

} 

ps[i] = make_pair(sw, sv); 

} 

5 

/ 去 除 多 余 的 元 素 

sort(ps, ps + (1 << n2)); 

intm = 1; 

for (int i= 1; i < 1 << n2; i++) { 
6 

if (ps[m - 1].second < ps[i].second) { 
ps[m++] = psli]; 

} 

} 

7 

/ 枚 举 后 半 部 分 并 求解 


ll res = 0; 


for (int i = 0; i < 1 << (n - n2); i++) { 
ll sw = 0, sv = 0; 

for (int j = 0; j <n - n2; j++) { 
8 

if (i>>j &1){ 

sw += w[n2 + j]; 

sv += v[n2 + j]; 

} 

} 

if (sw <= W) { 

9 


ll tv = (lower_bound(ps, ps + m, make_pair(W - sw, INF)) - 1)->second; res 
= max(res, sv + tv); 


} 

} 

10 

printf("%lld\n", res); 
} 

11 


i 1 
5 


P 


AT/ ++ I 
A DO 


164% 3 Æ 


3.2.5 坐标 离散 化 
区 域 的 个 数 


wx h 的 格子 上 画 了 n 条 或 牌 直 或 水 乎 的 宽度 为 1 的 直线 。 求 出 这 些 线 
将 格子 划分 成 了 多 少 个 区 域 。 


例 
限制 条 件 


[O 1 <w, h < 1000000 


[}1<n<500 
样 例 


w=10,h=10,n=5 


xl = {1, 1, 4, 9, 10} 
x2 = {6, 10, 4, 9, 10} 
yl = {4, 8, 1, 1,6} 

y2 = {4, 8, 10, 5, 10} 
(对 应 于 前 面 的 例 图 ) 


准备 好 wxh 的 数组 ， 并 记录 是否 有 直线 通过 ， 然 后 参考 2.1 市 利用 深 
度 优 先 搜索 求 水 寺 数 的 方法 ， 


可 以 求 出 被 分 割 出 的 区 域 个 数 。 但 是 这 个 问题 中 w 和 hh 最 大 为 
1000000， 所 以 没 办 法 创建 w x h 的 数 


组 。 因 此 我 们 要 使 用 坐标 离散 化 这 一 技巧 。 

坐标 离散 化 示例 

3.2 常用 技巧 精 选 (一 ) 165 

如 上 图 所 示 ， 将 前 后 没有 变化 的 行列 消除 后 并 不 会 影响 区 域 的 个 数 。 
1 


数组 里 只 需要 存储 有 直线 的 行列 以 及 其 前 后 的 行列 就 足够 了 ， 这 样 的 
话 大 小 最 多 6 n x6n WET ° 


因此 就 可 以 创建 出 数组 并 利用 搜索 山 求 出 区 域 的 个 数 了 。 
2 


/输入 


int W, H, N; 
int X1[MAX_N], X2[MAX_N], YI[MAX_N], Y2[MAX_N |; 
/ 填充 用 


bool fld[MAX_N * 6][MAX_N * 6]; 

3 

/对 x1 和 xX2 进 行 坐标 离散 化 ， 并 返回 离散 化 之 后 的 视 度 
int compress(int *x1, int *x2, int w) { 

vector<int> xs; 

4 

for (int i = 0; i < N; i++) { 

for (int d = -1; d <= 1; d++) { 

int tx1 = x1[i] + d, tx2 = x2[i] + d; 


if (1 <= tx1 && txl <= W) xs.push_back(tx1); if (1 <= tx2 && tx2 <= W) 
xs.push_back(tx2); 5 


} 

} 

sort(xs.begin(), xs.end()); 
xs.erase(unique(xs.begin(), xs.end()), xs.end()); 

6 

for (int i = 0; i < N; i++) { 

x1[i] = find(xs.begin(), xs.end(), x1[i]) - xs.begin(); 
x2[i] = find(xs.begin(), xs.end(), x2[i]) - xs.begin(); 
} 

7 


return xs.size(); 

} 

void solve() { 

8 

/ 坐标 离散 化 

W = compress(X1, X2, W); 

H = compress(Y1, Y2, H); 
/填充 有 直线 的 部 分 

memset(fld, 0, sizeof(fld)); 

9 

for (int i = 0; i < N; i++) { 

for (int y = Y1[i]; y <= Y2[i]; y++) { 
for (int x = X1[i]; x <= X2[i]; x++) { 
fld[y][x] = true; 

} 


10 


11 


O 区 域 可 能 很 大 ， 所 以 用 递归 函数 实现 的 话 可 能 会 栈 洲 出 。 
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1/ 求 区 域 的 个 数 

int ans = 0; 

for (int y = 0; y < H; y++) { 

for (int x = 0; x < W; x++) { 

if (fldly][x]) continue; 

ans++; 

/ 宽度 优先 搜索 

queue<pair<int, int> > que; 

que.push(make_pair(x, y)); 

while (!que.empty()) { 

int sx = que.front().first, sy = que.front().second; 

que.pop(); 

for (int i = 0; i < 4; i++) { 

int tx = sx + dx[i], ty = sy + dylil; 

if (tx < 0 || W <= tx || ty < 0 || H <= ty) continue; if (fld[ty][tx]) continue; 
que.push(make_pair(tx, ty)); 

fld[ty][tx] = true; 

} 


intf("%d\n", ans); 
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4.3 成 为 图 论 大 师 之 路 


图 是 非常 有 用 的 数据 结构 ， 除 了 在 第 2 章 已 经 介绍 的 算法 外 ， 还 有 各 种 
各 样 的 相关 算法 。 在 此 ， 


我 们 主要 讨论 强 连通 分 量 分 解 和 最 近 公 共 人 祖先 等 问题 。 
4.3.1 强 连通 分 量 分 解 


对 于 一 个 有 回 图 顶点 的 子 集 S ， 如 末 在 5 内 任 取 两 个 顶点 U 和 vv ， 都 
能 找到 一 条 从 u 到 v 的 路 径 ， 那 么 


就 称 S 是 强 连 通 的 。 如 果 在 强 连 通 的 顶点 集合 S 中 加 入 其 他 任意 顶点 
集合 后 ， 它 部 不 再 是 强 连通 的 ， 


那么 就 称 S 是 原 图 的 一 个 强 连通 分 量 (SCC: Strongly Connected 
Component) 。 任 意 有 问 图 都 可 以 


分 解 成 看 干 不 相交 的 强 连通 分 量 ， 这 融 是 强 连 通 分 量 分 解 。 把 分 解 后 
的 强 连 通 分 量 缩 成 一 个 顶点 ， 


就 得 到 了 一 个 DAG (有 向 无 环 图 ) 
虚线 包围 的 部 分 构成 一 个 强 连通 分 量 


强 连通 分 量 分 解 可 以 通过 两 次 简单 的 DFS 实 现 。 第 一 次 DFS 时 ， 远 取 任 
意 顶点 作为 起 点 ， 授 历 所 


有 尚未 访问 过 的 顶点 ， 并 在 回调 前 给 顶点 标号 (post order, Jari 
历 ) 。 对 剩余 的 未 访问 过 的 顶 


点 ， 不 断 重 复 上 述 过 程 。 


完成 标号 后 ， 越 接近 图 的 尾部 (搜索 树 的 叶子 )  ， 顶 点 的 标号 越 小 。 
第 二 次 DFS 时 ， 先 将 所 有 边 


反问 ， 然 后 以 标号 最 大 的 顶点 为 起 点 进行 DFS。 这 样 DFS 所 通 历 的 顶点 
集合 惑 构成 了 一 个 强 连 通 


分 量 。 之 后 ， 只 要 还 有 疝 未 访问 的 顶点 ， 了 束 从 中 选取 标号 最 大 的 顶 辟 
不 断 重 复 上 述 过 程 。 


MN. 
N. 
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4.3 成 为 图 论 大 师 之 路 321 

1 

2 

后 续 遍 历 的 例子 。 根 据 搜索 顺序 的 不 同 ， 标 号 结果 也 可 能 不 同 
3 


4 
反问 后 的 图 


正如 前 文 所 述 ， 我 们 可 以 将 强 连 通 分 量 缩 点 并 得 到 DAG。 此 时 可 以 发 
现 ， 标 号 最 大 的 万 点 焉 属于 


5 


DAG 头 部 (搜索 树 的 根 ) MEN E ° A, RL AA, HARE 
沿边 访问 到 这 个 强 连通 分 量 


以 外 的 顶点 。 而 对 于 强 连 通 分 量 内 的 其 他 顶点 ， 其 可 达 性 不 受 边 反问 
的 影响 ， 因 此 在 第 二 次 DFS 


时 ， 我 们 可 以 遍历 一 个 强 连 通 分 量 里 的 所 有 项 点。 

6 

7 

8 

9 

边 反 向 后 ， 从 8、9、10 号 顶点 只 能 到 达 其 头 部 方向 的 顶点 11 和 12 
该 算法 只 进行 了 两 次 DFS， 因 而 总 的 复杂 度 是 O (VEE) > 
int V; // 顶点 数 

10 


vector<int> G[MAX_V]; // 图 的 邻接 表 表 示 


vector<int> rG[IMAX_V]; // 把 边 反 向 后 的 图 
vector<int> vs; // 后 序 涡 历 顺 序 的 顶点 列表 
bool used[MAX_V]; // 访问 标记 

int cmp[MAX_V]; / 所 属 强 连通 分 量 的 拓扑 友 


11 
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void add_edge(int from, int to) { 
G[from].push_back(to); 
rG[to].push_back(from); 

} 

void dfs(int v) { 

used[v] = true; 

for (int i = 0; i < G[v].sizeQ; i++) { 
if (tused[G[v][i]]) dfs(G[vI[iD; 

} 

vs.push_back(v); 

} 

void rdfs(int v, int k) { 

used[v] = true; 

cmplv] = k; 

for (int i = 0; i < rG[v].sizeQ; i++) { 
if (tused[rG[v][i]]) rdfs(rG[v][il, k); 
} 

} 


int scc() { 


memset(used, 0, sizeof(used)); 
vs.clear(); 

for (int v = 0; v < V; v++) { 

if (!used[v]) dfs(v); 

} 

memset(used, 0, sizeof(used)); 
int k = 0; 

for (int i = vs.size() - 1; i >= 0; i--) { 
if (tused[vs[i]]) rdfs(vs[i], k++); 
} 

return k; 

j 

Popular Cows (POJ No.2186) 


每 头 牛 都 想 成 为 牛 群 中 的 红 人 “。 给 定 N 头 牛 的 牛 群 和 M 个 有 序 对 ( A， 
B)。(A,B) 表 示 牛 A 认 


为 牛 B 是 红 人 。 该 关系 具有 传递 性 ， 所 以 如 果 牛 A 认为 牛 B 是 红 人 ， 
牛 B 认为 牛 C 是 红 人 ， 


那么 牛 A 也 认为 牛 C 是 红 人 。 不 过 ， 给 定 的 有 序 对 中 可 能 包含 (A,B) 
和 ( B,C)， 但 不 包含 (A,C)。 


求 被 其 他 所 有 牛 认为 是 红 人 的 牛 的 总 数 。 
限制 条 件 


[] 1 <N < 10000 


[] 1 < M < 50000 
[1<A,B<N 

4.3 成 为 图 论 大 师 之 路 323 
样 例 

1 

输入 

N=3 

M=3 

(A, B) = {(1, 2), (2, 1), 2, 3)} 
2 

输出 


著 虚 以 牛 为 顶点 的 有 向 图 ， 对 每 个 有 序 对 ( A , B ) 连 一 条 从 A 到 B 的 有 
回 边 。 那 么 ， 被 其 他 所 有 和 牛 认 


为 是 红 人 的 牛 对 应 的 顶 扎 ， 也 束 是 从 其 他 所 有 顶点 都 可 达 的 顶点 。 虽 
然 这 可 以 通过 从 每 个 顶点 出 


发 搜索 求 得 ， 但 总 的 复杂 度 却 是 O (NM )， 是 不 可 行 的 ， 必 须要 考虑 
更 为 高 效 的 算法 。 


4 


假设 有 两 头 牛 A 和 B 都 被 其 他 所 有 和 牛 认为 是 红 人 “。 那 么 显然 ，A 被 B 认 为 
ELA, BUCAY NEL 


人 ， 即 存在 一 个 包含 A、B 两 个 顶点 的 圈 ， 或 者 说 ，A、B 同 属于 一 个 强 
连通 分 量 。 反 之 ， 如 果 一 

aah 那么 其 所 属 的 强 连通 分 量 内 的 所 有 人 牛 
ARB AT PUA E 

5 

人 。 由 此 ， 我 们 把 图 进行 强 连 通 分 量 分 解 后 ， 至 多 有 一 个 强 连通 分 量 


满足 题目 的 条 件 。 而 按 前 面 


介绍 的 算法 进行 强 连通 分 量 分 解 时 ， 我 们 还 能 够 得 到 各 个 强 连通 分 量 


拓扑 排序 后 的 顺序 ， 唯一 可 
能 成 为 解 的 只 有 拓扑 序 最 后 的 强 连通 分 量 。 所 以 在 最 后 ， 我 们 只 要 检 


得 这 个 强 连通 分 量 是 否 从 所 


6 


有 顶点 可 达 就 好 了 。 该 算法 的 复杂 度 为 O(N+ M)， 足 以 在 时 限 内 解 
决 原 题 。 


/输入 


int N, M; 


int AAMAX_M], BIMAX_M]; 


7 


void solve() { 


V=N; 


for (int i = 0; i < M; i++) { 


add_edge(A[i] - 1, Bli] - 1); 


8 


} 
int n = scc(); 
统计 备 选 解 的 个 数 
int u = 0, num = 0; 
9 
for (int v = 0; v < V; v++) { 
if (cmp[v] == n - 1) { 
U= y; 
num++; 
} 
10 
} 
1/ 检查 是 否 从 所 有 后 可 达 
memset(used, 0, sizeof(used)): 
rdfs(u, 0); // 重用 强 连 通 分 量 分 解 的 代码 


11 
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for (int v = 0; v < V; v++) { 
if (!used[v]) { 
/ 从 该 点 不 可 达 


printf("%d\n", num); 
} 
4.3.2 2-SAT 


给 定 一 个 布尔 方程 ， 判 断 是 否 存在 一 组 布尔 变量 的 真 值 指派 使 整个 方 
程 为 真 的 问题 ， 被 称 为 布 


尔 方程 的 可 满足 性 问题 (SAT) 。SAT 问 题 是 NP 完全 的 ， 但 对 于 满足 一 
定 限制 条 件 的 SAT 问 题 ， 还 


是 能 够 有 效 求解 的 。 我 们 将 下 面 这 种 布尔 方程 称 为 合 取 范 式 。 
(aVDV...)A(CV dV...)A... 


其 中 a,b ,.…. 称 为 文字 ， 它 十 一 个 布尔 变量 或 其 否定 。 像 (a VD V.…) 这 
样 用 Vv 连接 的 部 分 称 为 子 句 。 


如 琳 合 取 范 式 的 每 个 子 句 中 的 文字 个 数 都 不 超过 两 个 ， 那 么 对 应 的 SAT 
问题 又 称 为 2-SAT 问 题 。 


O 2-SAI 布 尔 公 式 的 例子 
D(aVb)A-~a 令 a 为 假 而 b 为 真 ， 则 可 以 满足 


D(aV-b)A(bVc)A(-cV-a) 令 a 和 4b 为 真 而 c 为 假 ， 则 可 以 满 
E 


D(aVb)A(aV-b)A (maVb)A (~-avV-b) 无 法 满足 


利用 强 连通 分 量 分 解 ， 可 以 在 布尔 公式 子 句 数 的 线性 时 间 内 解决 2-SAT 
问题 。 首 先 ， 利 用 = (A 


W) 将 每 个 子 句 ( a Vb ) 改 写成 等 价 形式 (J a => b mb > a )°。 这 样 原 布 
尔 公 式 就 变 成 了 把 a => b 形式 的 


布尔 公式 用 人 连接 起 来 的 形式 。 对 每 个 布尔 变量 x ， 构 造 两 个 顶点 分 别 
代表 x 和 - x ， 以 = 天 系 为 边 


建立 有 向 图 。 此 时 ， 如 果 图 上 的 a 点 能 够 到 达 b 点 的 话 ， 就 表示 当 a 
为 真 时 也 一 定 为 真 。 因 此 ， 该 


图 中 同一 个 强 连通 分 量 中 所 偏 的 所 有 文字 的 布尔 值 均 相同 。 


如 果 存 在 某 个 布尔 变量 x ，x Alo x 均 在 同一 个 强 连通 分 量 中 ， 则 显然 
无 法 令 整 个 布尔 公式 的 值 为 真 。 


有 反之， 如 末 不 存在 这 样 的 布尔 变量 ， 那 么 对 于 每 个 布尔 变量 x ， 让 
x 所 在 的 强 连通 分 量 的 拓扑 序 在 -x 所 在 的 强 连通 分 量 之 后 x A 
号 是 使 得 该 公式 的 值 为 真 的 一 组 合适 的 布尔 变量 赋值 。 


(av-b)AbvaA(-cv=a) 所 对 应 的 图 

int main() { 

/布尔 公式 为 (av =b)A (DV OA (acV -al 时 
/ 构造 6 个 顶点 ， 分 别 对 应 a\`b、c、-a、 ab > ace e 
4 

V =6; 

/ av bé i-a=>-bAb=a 

add_edge(3, 4); // 从 -na 连 一 条 到 -的 边 
add_edge(1, 0); // 从 b 连 一 条 到 a 的 边 

5 

/byv c 转 成 "b=:cA-mc=b 

add_edge(4, 2); / 从 -b 连 一 条 到 c 的 边 
add_edge(5, 1); // 从 -~c 连 一 条 到 b 的 边 

/ -CV na 转 成 cnaAa 一 -C 

add_edge(2, 3); / 从 c 连 一 条 到 -a 的 边 

6 

add_edge(0, 5); / 从 a 连 一 条 到 -~c 的 边 


/ 进行 强 连通 分 量 分 解 


scc(); 


/ 判断 十 否 x 和 -x 在 不 同 的 强 连通 分 量 中 


7 

for (int i = 0; i < 3; i++) { 
if (cmp[i] == cmp[3 + i]) { 
printf(""NO"); 

return 0; 

} 

8 

} 

/ 如 有 宁可 满足 ， 则 给 出 一 组 解 
printf("YES\n"); 

for (int i = 0; i < 3; i++) { 
9 

if (cmp[i] > cmp[3 + i]) { 
printf("true\n"); 

} else { 

printf("false\n"); 

} 

10 

} 


return 0; 


11 


高 级 篇 
Priest John’s Busiest Day (POJ No.3683) 


约翰 是 街区 里 唯一 的 神父 。 假 设 有 NN 对 新 人 打算 在 同一 天 举行 结婚 仪 
式 。 第 i 对 新 人 的 结婚 


仪式 的 时 间 为 5i 到 Ti ， 在 其 仪式 开始 时 或 是 结束 时 需要 进行 一 个 用 时 
Di 的 特别 仪式 (也 就 


是 从 Si 到 Si+ Di 或 是 从 Ti0 Di 到 Ti) ， 该 特别 仪式 需要 神父 在 场 。 
请 判断 是 否 可 以 通过 合理 安 


排 每 个 特别 仪式 在 开始 还 是 结束 时 进行 ， 从 而 保证 神父 能 够 出 席 所 有 
的 符 别 仪式 。 如 果 可 能 


的 话 ， 请 输出 出 席 各 个 特别 仪式 的 时 间 。 当 然 ， 神 父 不 可 能 同时 出 席 
多 个 特别 仪式 。 不 过 神 


父 前 往 仪 式 的 途中 所 花费 的 时 间 可 以 名 略 不 计 ， 神 父 可 以 在 出 席 完 一 
个 特别 仪 后 ， 立 刻 出 席 


另 一 个 开始 时 间 与 其 结束 时 间 相 等 的 特别 仪式 。 
限制 条 件 


[] 1< N <1000 
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样 例 

输入 

N=2 

(S, T, D) = {(08:00, 09:00, 30), (08:15, 09:00, 20)} 
输出 

YES 

08:00 08:30 

08:40 09:00 


对 于 每 个 结婚 仪式 i ， 只 有 在 开始 或 结束 时 进行 特别 仪式 两 种 选择 。 
此 可 以 定义 变量 xi 


xi 为 真 e 在 开始 时 进行 特别 仪式 


这 样 ， 对 于 结婚 仪式 i 和 j， 如 果 Si~ Si+ Di 和 Sj~ Sj+ Di 冲突， 就 有 
-xi Va xj AR CN FIIRT R ` a 


束 和 开始 、 结 束 和 结束 等 三 种 情况 ， 也 可 以 得 到 类 似 的 条 件 。 于 是 ， 
要 保证 所 有 特别 仪式 的 时 间 


不 冲突 ， 只 要 考虑 将 这 所 有 的 子 句 用 人 连接 起 来 所 得 到 的 布尔 公式 整 好 
了 。 例 如 ， 对 于 输入 样 例 ， 


可 以 的 到 布尔 公式 


(-x1Vv-x2)A(x1V-x2)A(x1VX2) 而 当 x1 为 真 而 x 2 为 假 时 ， 其 值 
为 真 。 这 样 ， 我 们 就 把 原 问 题 转 为 了 2-SAT 问 题 。 接 下 来 只 要 进行 强 


使 得 布尔 公式 值 为 真 的 一 组 布尔 变量 峰值 
RLF T° 


/输入 


int N; 

int SIMAX_N], TIMAX_N], DIMAX_N]; // S 和 T 是 换算 成 分 钟 后 的 时 间 
4.3 成 为 图 论 大 师 之 路 327 

1 

void solve() { 

//Q~N-1: xi 

// N2N-1: =x_1 

V=N *2; 

for (int i = 0; i < N; i++) { 

for (int j = 0; j < i; j++) { 

2 

if (min(S[i] + Dil, S[j] + DIJI) > max(S[il, SID) { 
// Xi 一 -Xj]、X j=7X 1 

add_edge(i, N + j); 

add_edge(j, N + i); 

} 

3 

if (min(S[i] + Di], THI) > max(S[i], TH] - DGD) { 
// X 1>X_j ` XL 


add_edge(i, j); 


add_edge(N +j, N + i); 

} 

4 

if (min(T[i], S[j] + DII) > max(T[i] - D[il, SID) { 
// ~X_ i=7X ` X_J=x_i 

add_edge(N + i, N + j); 

add_edge(j, i); 

} 

5 

if (min(T[il, TGD) > max(T[i] - D[il, TE] - DD) { 
// 7X_1=X_} ` -x jx i 

add_edge(N + i, j); 

add_edge(N + j, i); 

} 

6 

} 

} 

scc(); 

/判断 是 否 可 满足 


7 


for (int i = 0; i < N; i++) { 

if (cmp[i] == cmp[N + iJ) { 
printf("NO\n"); 

return; 

} 

} 

8 

/ 如 果 可 满足 ， 则 给 出 一 组 解 
printf("YES\n"); 

for (int i = 0; i < N; i++) { 

if (cmp[i] > cmp[N + i]) { 

9 

/x_i 为 真 ， 即 在 结婚 仪式 开始 时 举行 


printf("%02d:%02d %02d:%02d\n", S[i] / 60, S[i] % 60, (S[i] + D[i]) / 60, 
(S[i] + D[i]) % 60); 


} else { 
// x_i 为 假 ， 即 在 结婚 仪式 结束 时 举行 
10 


printf("%02d:%02d %02d:%02d\n", (T[i] - D[i]) / 60, (T[i] - DID) % 60, 
T[i] / 60, T[i] % 60); 


} 


11 
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4.3.3 LCA 


在 有 根 树 中 ， 两 个 节点 u 和 v 的 公共 祖先 中 距离 最 近 的 那个 被 称 为 最 
近 公 共 祖 先 (LCA, Lowest Common Ancestor)。 用 于 高 效 计算 LCA 的 
算法 有 许多 ， 在 此 我 们 介绍 其 中 的 两 种 。 在 下 文中 ， 我 

们 都 假设 节点 数 为 n。 

LCA 的 例子 (4 和 7 的 LCA 为 2，8 和 6 的 LCA 为 1，5 和 8 的 LCA 为 5) 

1. 基于 二 分 搜索 的 算法 


记 节 点 v 到 根 的 深度 为 depth(v)。 那 么 ， 如 果 节 点 w 是 u 和 Vv 
的 公共 祖先 的 话 ， 让 u 向 上 走 


depth( u )-depth( w ) 步 ， 让 vy 同上 走 depth(v )-depth( w ) 步 ， 束 都 将 走 到 
w。 因 此， 首先 让 wu 和 vy 中 较 深 的 


一 方 同上 走 |depth( u )-depth(v )| 步 ， 再 一 起 一 步 步 辐 上 走 ， 直 到 走 到 同 
一 个 万 点 ， 束 可 以 在 


O (depth( u )+depth(v )) 时 间 内 求 出 LCA 。 
/ 输入 

vector<int> GIMAX_V]; // 图 的 邻接 表 表 示 
int root; // 根 斑点 的 编号 


int parent[IMAX_V]; // 父亲 节点 〈 根 节点 的 父亲 记 为 -1) 


int depth[MAX_V]; // 节点 的 深度 
void dfs(int v, int p, int d) { 
parent[v] = p; 

depth[v] = d; 

for (int i = 0; i < G[v].sizeQ; i++) { 


if (G[v]Li] != p) dfs(G[v]Lil, v, d + 1); 


} 

} 

/ 预 处 理 

void init() { 

// 预 处 理 出 parent 和 depth 

dfs(root, -1, 0); 

} 

/ 计算 u 和 v 的 LCA 
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int Ica(int u, int v) { 

1 

// 让 uv 向 上 走 到 同一 深度 

while (depth[u] > depth[v]) u = parent[u]; 
while (depth[v] > depth[u]) v = parent[v]; 
avia LER — A 

while (u != v) { 

u = parent[u]; 

2 

v = parent[v]; 


} 


return U; 
} 
3 


节点 的 最 大 深度 是 O(n)， 所 以 该 算法 的 复杂 度 也 是 O(n )。 如 果 只 需 
计算 一 次 LCA 的 话 ， 这 便 足 够 


了 。 但 如 果 要 计算 多 对 节点 的 LCA 的 话 如 何 是 好 呢 ? 刚才 的 算法 ， 通 
LAA LR ARA 


计算 wu 和 vv 的 LCA。 这 里 ， 到 达 了 同一 节点 后 ， 不 论 再 怎么 同上 走 ， 
PITA Si Aide I] A All 


4 


用 这 一 点 ， 我 们 能 够 利用 二 分 搜索 求 出 到 达 共 同 祖先 所 需 的 最 少 步 数 
吗 ? 事实 上 ， 只 要 利用 如 下 


预 处 理 ， 束 可 以 实现 二 分 搜索 。 


首先 ， 对 于 任意 顶点 v ， 利 用 其 父 杀 节点 信息 ， 可 以 通过 parent2[ v 
]=parent[parent[ v 得 到 其 同上 走 


5 


两 步 所 到 的 顶点 。 再 利用 这 一 信息 ， 又 可 以 通过 parent4[ v 
]=parent2[parent2[ v ]] 得 到 其 问 上 走 四 步 


所 到 的 顶点 。 依 此 类 推 ， 残 能 够 得 到 其 问 上 走 2K 步 所 到 的 顶点 parent[ 
k]iv]° ÆT k=floor(log mn) 以 内 


的 所 有 信息 后 ， 承 可 以 二 分 搜索 了 ， 每 次 的 复杂 上 度 是 (logn)? 5 
外 ， 预 处 理 parent[K][v] 的 复杂 度 


6 


是 O(nlogn)。 


/输入 

vector<int> G[MAX_V]; // 图 的 邻接 表 表 示 
int root; // AR PD AA Fa 

7 


int parent[MAX_LOG_V][MAX_V]; // 向 上 走 2^k 步 所 到 的 节点 (超过 根 
时 记 为 -1) 


int depth[MAX_V]; // 节点 的 深度 
void dfs(int v, int p, int d) { 

8 

parent[0][v] = p; 

depth[v] = d; 

for (int i = 0; i < G[v].sizeQ; i++) { 
if (G[v]Li] != p) dfs(G[v][il, v, d + 1); 
} 

} 

9 

/ 预 处 理 

void init(int V) { 

// 预 处 理 出 parent[0] 和 depth 
dfs(root, -1, 0); 


10 


// 预 处 理 出 parent 

for (int k = 0; k+ 1 < MAX_LOG_V; k++) { 
for (int v = 0; v < V; v++) { 

if (parent[k][v] < 0) parent[k + 1][v] = -1; 

else parent[k + 1][v] = parent[k][parent[k][v]]; 
11 


} 


AUCH 


加 
[CO 


/ 计算 u 和 v 的 LCA 


int Ica(int u, int v) { 

// 让 uv 向 上 走 到 同一 深度 

if (depth[u] > depth[v]) swap(u, v); 

for (int k = 0; k < MAX_LOG_V; k++) { 
if ((depth[v] - depth[u]) >> k & 1) { 

v = parent[k][v]; 

} 

} 

if (u == v) return u; 


/利用 二 分 搜索 计算 LCA 


for (int k = MAX_LOG_V - 1; k >= 0; k--) { 
if (parent[k][u] != parent[k][v]) { 

u = parent[k][u]; 

v = parent[k][v]; 

} 

} 

return parent[O][u]; 

} 


像 这 样 ， 预 处 理 出 2K 的 表 的 技巧 ， 在 计算 LCA 之 外 也 很 有 用 ， 相 天 的 
问题 也 经 常 出 现在 程序 设计 


竞赛 当中 。 对 此 ， 下 一 市 中 还 会 介绍 其 他 例子 。 
2. 基于 RMQ 的 算法 


对 于 涉及 有 根 树 的 问题 ， 将 树 转 为 从 根 DFS 标 号 后 得 到 的 序列 处 理 的 技 
UVF te TTARONT 


LCA， 利 用 该 技巧 也 能 够 高 效 地 计算 。 首 移 ， 按 从 根 DFS 访 问 的 顺序 得 
到 顶点 序列 vs[ 让 和 对 应 的 


深度 depth[ i]。 对 于 每 个 顶点 v ， 记 其 在 vs 中 首次 出 现 的 下 标 为 id[ v 
] o 


样 例 对 应 的 标号 


这 些 都 可 以 在 O(n ) 时 间 内 求 得 。 而 LCA( u,v ) 束 是 访问 u 之 后 到 访问 
v 之 前 所 经 过 顶点 中 离 根 最 近 的 


那个 ， 假 设 id[ u ]<id[ v]， 那 么 有 


LCA (u,v)=vs[id[u]<i<id[v] 中 令 depth (i 最 小 的 门 


而 这 可 以 利用 RMQ 高 效 地 求 得 。 


4.3 成 为 图 论 大 师 之 路 331 
/ 输入 
1 


vector<int> G[MAX_V]; // 图 的 邻接 表 表 示 
int root; 

int vs[MAX_V * 2 - 1]; // DFS 访 问 的 顺序 
int depth[MAX_V * 2 - 1];/ 节点 的 深度 

2 

int id[IMAX_V];V 各 个 顶点 在 vs 中 首次 出 现 的 下 标 
void dfs(int v, int p, int d, int &k) { 

id[v] = k; 

vs[k] = v; 

3 

depth[k++] = d; 


for (int i = 0; i < G[v].size(); i++) { 


if (G[vJ[] != p) { 
dfs(G[v][i], v, d + 1, k); 
vs[k] = v; 

depth[k++] = d; 

4 

} 

} 

} 

// 预 处 理 

5 

void init(int V) { 

/ 预 处 理 出 vs、depth 和 id 
int k = 0; 

dfs(root, -1, 0, k); 

/ 预 处 理 出 RMQ (返回 的 不 是 最 小 值 ， 而 是 最 小 值 对 应 的 下 标 ) 
6 

rmq_init(depth, V * 2 - 1); 
} 

/ 计算 u 和 v 的 LCA 


int Ica(int u, int v) { 


7 

return vs[query(min(id[u], id[v]), max(id[u], id[v]) + 1)]; 

} 

8 

Housewife Wind (POJ No.2763) 

xx 村 里 有 nE, 7) eZ TA CT) GAARA, Pr) EN Ae 
一 棵 树 。 通 过 连接 ai 号 小 屋 和 bi 号 小 屋 的 道路 i 需要 花费 wi 的 时 间 。 
你 一 开始 在 s 号 小 屋 。 请 处 理 以 下 gq 个 查询 。 

9 


A: 输出 从 当前 位 置 移动 到 市 点 x 所 需 的 时 间 。B: 将 通过 道路 x 所 需 
的 时 间 改 为 +。 


限制 条 件 


[] 1< n <100000 


10 

[] 0< q <100000 
[]1<ai,bi<n 
[] 1< wi <10000 


11 


RAR 
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样 例 
输入 


n=3 


q=3 
s=1 

(a, b, w) = {(1, 2, 1), (2, 3, 2)} 

(查询 的 类 型 , x(, t)) = {(A, 2), (B, 2, 3), (A, 3)} 
输出 

从 1 移动 到 2 

从 2 移动 到 3 


虽然 直接 DFS 也 可 以 求 出 树 上 两 点 之 间 的 距离 ， 但 是 这 对 于 每 个 A 类 型 
的 查询 ， 都 要 花费 O(n ) 的 


时 间 ， 实 在 太 慢 了 。 必 须 利 用 树 的 特性 ， 得 到 更 为 高 效 的 算法 。 为 了 
高 效 地 处 理 A 类 型 的 查询 ， 


可 以 利用 二 分 搜索 版 LCA 算 法 中 介绍 的 技巧 ， 记 杂 下 从 每 个 顶 后 同上 
走 2K 步 的 总 长 度 。 这 样 一 来 ， 


在 O (dogn) 地 计算 LCA 的 同时 ， 也 可 以 同样 O (log n ) 地 求 出 到 LCA 的 
距离 ， 因 此 处 理 A 类 型 查询 的 


复杂 度 为 (log n )。 但 是 ， 这 个 方法 对 于 B 类 型 的 查询 却 无 法 高 效 地 处 
理 。 


因此 ， 让 我 们 先 考虑 一 下 图 是 链 状 时 这 一 简单 的 情况 。 假 设 i 和 i+1 两 
扩 之 间 的 边 的 长 度 为 wi ， 则 


两 点 u 和 v(u<v) 之 间 的 距离 为 
v1 

O 

DO wi 

¡Du 


只 要 用 BIT， 不 论 是 A 类 型 的 查询 还 是 B 类 型 的 查询 ， 都 能 够 在 (log n 
) 时 间 内 处 理 > 


在 回 过 头 来 考虑 树 的 情况 。 因 为 树 中 连接 两 点 的 路 径 是 唯一 的 ， 如 果 
我 们 对 顶点 进行 合理 排列 的 


话 ， 能 人 否 像 链 状 时 那样 ， 进 行 类 似 的 处 理 呢 ? 考虑 利用 RMQ 计 算 LCA 
时 所 用 的 ， 按 DFS 访 问 的 顺 


序 排列 的 顶点 序列 。 这 样 ，u 和 v 之 间 的 路 径 ， 就 是 在 序列 中 u 和 v 
之 间 的 所 有 边 减 去 往返 重复 的 部 


分 得 到 的 结果 。 
将 树 转 为 序列 
4.3 成 为 图 论 大 师 之 路 333 


于 是 ， 只 要 令 边 的 权重 沿 叶 子 方 同 为 正 ， 沿 根 方 同 为 负 ， 那 么 往返 重 
复 的 部 分 就 目 然 抵消 了 ， 于 


1 
是 有 


(u,v 之 间 的 距离 )=( 从 LCA (u,v ) 到 4w 的 边 的 权重 和 )+( 从 LCA (u,v) 
到 v 的 边 的 权重 和 ) 2 


同 链 状 的 情况 一 样 ， 利 用 BIT 的 话 ， 计 算 权 重 和 和 更 新 边 权 都 可 以 在 O 
(ogn ) 时 间 内 办 到 ， 而 LCA 也 能 够 在 O (log n ) 时 间 内 求 得 。 


struct edge { int id, to, cost; }; 

3 

int n, q, S; 

int fMAX_V - 1], bIMAX_V - 1], wWIMAX_V - 1]; 
int type[MAX_Q]; //0: A 类 型 ，1: BR 
intx[MAX_Q], t[MAX_Q]; 

4 

vector<edge> G[MAX_V]; // 图 的 邻接 表 表 示 

int root; 

int Vs[MAX_V * 2 - 1]; // DFS 访 问 的 顺序 


int depth[MAX_V * 2 - 1]; / 节点 的 深度 


5 
int idIMAX_V];V/ 各 个 顶点 在 vs 中 首次 出 现 的 下 标 


int es[((MAX_V - 1) * 2]; / 边 的 下 标 (i*2+( 叶 子 方向 :0, 根 方向 :1)) 


void dfs(int v, int p, int d, int &k) { 
id[v] = k; 

6 

vs[k] = v; 

depth[k++] = d; 

for (int i = 0; i < G[v].size(); i++) { 


edge «ze = G[v]L[i]; 


if (e.to != p) { 
7 
add(k, e.cost); 


es[e.id * 2] = k; 
dfs(e.to, v, d + 1, k); 
vs[k] = v; 
depth[k++] = d; 
add(k, -e.cost); 

8 


es[e.id * 2 + 1] = k; 


9 

int stack_v[MAX_V + 10]; 

int stack_i[MAX_V + 10]; 

/ 预 处 理 

void init(int V) { 

10 

/ 初始 化 BIT 

bit_n = (V - 1) * 2; 

// 预 处 理 出 vs、depth、id 和 es 
int k = 0; 

dfs(root, -1, 0, k); 

11 

/ 预 处 理 出 RMQ (返回 的 不 是 最 小 值 ， 而 是 最 小 值 对 应 的 下 标 ) 
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rmq_init(depth, V * 2 - 1); 
} 
/ 计算 u 和 v 的 LCA 


int Ica(int u, int v) { 

return vs[query(min(id[u], id[v]), max(id[ul, id[v]) + 1)]; 
} 

void solve() { 

/ 预 处 理 

root = n/ 2; // 不 论 以 哪个 节点 为 根部 没有 问题 

for (int i= 0;i<n-1;it+){ 

Glali] - 1].push_back((edge){i, bli] - 1, wlil}); 


G[b[i] - 1].push_back((edge){i, ali] - 1, wlil}); 


} 
init(n); 
/ 处 理 查 询 


intv=s-1;V/ 当 前 位 置 
for (int i = 0; i < q; i++) { 
if (typeli] == 0) { 

/ 从 当前 位 置 移动 到 xj 
int u = x[i] - 1; 

int p = Ica(v, u); 


/利用 BIT 计算 p 到 v 和 p 到 u 的 费用 之 和 ， 即 区 间 (id[p],id[v]] 和 
(id[pj,id[u]] 的 权重 和 


printf("%d\n", sum(id[v]) + sum(id[u]) - sum(id[p]) * 2); v = u; 


} else { 
HPS xX AEA Li] > 
int k = x[i] - 1; 

add(es[k * 2], t[i] - w[k]); 
add(es[k * 2 + 1], w[k] - t[i]); 
w[k] = tli]; 

) 

) 

) 


